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Abstract 

Null Polygon Wilson Loops inM — ^ SYM can be computed using the Operator Product 
Expansion in terms of a transition amplitude on top of a color Flux tube. That picture is 
valid at any value of the 't Hooft coupling and is studied here in the planar limit. So far it has 
been efficiently used at weak coupling in cases where only a single particle is flowing. At any 
finite value of the coupling however, an infinite number of particles are flowing on top of the 
color flux tube. A major open problem in this approach was how to deal with generic multi- 
particle states at weak coupling. In this paper we study the propagation of any number 
of flux tube excitations at weak coupling. We do this by flrst mapping the Wilson loop 
expectation value into a sum of two point functions of local operators. That map allows us to 
translate the integrability techniques developed for the spectrum problem back to the Wilson 
loop. In particular, we find that the fiux tube Hamiltonian can be represented as a simple 
kernel acting on the loop. Having an explicit representation for the flux tube Hamiltonian 
allows us to treat any number of particles on an equal footing. We use it to bootstrap 
some simple cases where two particles are flowing, dual to N^MHV amplitudes. The flux 
tube is integrable and therefore has other (inflnite set of) conserved charges. The generating 
function of all of these charges is constructed from the monodromy matrix between sides 
of the polygon. We compute it for some simple examples at leading order in perturbation 
theory. At strong coupling, these monodromies were the main ingredients of the Y-system 
solution. To connect the weak and strong coupling computations, we study a case where an 
infinite number of particles are propagating already at leading order in perturbation theory. 
We obtain a precise match between the weak and strong coupling monodromies. That match 
is the Wilson loop analog of the well known Prolov-Tseytlin hmit where the strong and weak 
coupling descriptions become identical. Hopefully, putting the weak and strong coupling 
descriptions on the same footing is the first step in understanding the all loop structure. 
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1 Introduction 
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Figure 1: The OPE picture. Flux tube excitations are created at the bottom and absorbed in the 
top. The flux is thick in spacetime but is 1+1 dimensional in AdS. Its excitations are integrable. 

The Operator Product Expansion (OPE) for Null Polygon Wilson Loops translates the 
Wilson loop computation into a transition amplitude on top of a color Flux tube [1], see 
figure 1. The OPE does not rely on perturbation theory. Indeed, it was observed both at 
weak and at strong coupling in TV = 4 SYM. The main hope of this program is to fill up the 
gap and find the same description at any intermediate coupling for this particular theory. 

The flux preserves some symmetries. We can therefore classify the states by the way they 
transform under these symmetries. In particular the states have energies and continuous 
momenta} The states are gapped at any value of the coupling [2,3] and therefore we can 
further classify them according to the number of particles. 

One main difference between weak and strong coupling is in the number of particles 
involved. At a given order in the weak coupling perturbative expansion, we have a maximum 
number of particles that can be excited and flow in the flux tube. Sometimes we can have 
at most a single excitation. This is the simplest case and is the one that is currently well 
understood from the OPE point of view.^ On the other hand, at strong coupling, we have a 
continuous distribution of infinitely many flux tube excitations. 

Hence, to further connect the weak and strong coupling we should go beyond the single 
^Also known as twist and conformal spin. 

^Despite the simplicity of the single particle exchange case, it turns out to be a very powerful tool for 
computing or constraining scattering amplitudes at the first few loop orders. For example, under simple 
assumptions, it fixes completely the two loops MHV amplitudes and the one loop NMHV amplitudes [4-6]. 
At higher loops it does not fix the full result but provides very helpful tight constraints on the result. 
They were recently used in constraining the hexagon NMHV amplitude at two loops [7], the hexagon MHV 
amplitude at three loops [8] and for the three loop 2D MHV octagon [9]. The single particle OPE was also 
used in deriving the Q equation [10], see also [11]. 
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Figure 2: At strong coupling the polygon Wilson loop expectation value is computed by 
a minimal surface area. The surface ends on the polygon at the boundary of AdS and is 
stretched in the AdS radial direction. The holonomy T{u) of the flat connection between two 
edges of the polygon is drawn in red. It is a measure of all the higher conserved charges of 
what is flowing through it. Such holonomies are the building blocks in the Y-system strong 
coupling solution [12]. 

particle OPE and get a better handle on the flux tube states. That is the main motivation 
for this paper. 

The reason why a single particle is considerably simpler is that single particle wave 
functions are always trivial to classify. They are just plane waves with definite momenta. On 
the other hand, many particles can interact. As a result, what the two or more particle wave 
functions are becomes a dynamical question that depends on the details of the interaction. 
To handle these states, we need to understand the map to the flux tube wave functions more 
concretely, as we will do in this paper. 

The flux tube states can be characterized by their charges. The most obvious of all is 
the energy. By reading the energy of these states we can easily make an infinite number 
of predictions at any loop order [1]. In sections 2-4 we focus on the Hamiltonian which 
measures this energy. We start by considering a simple set of examples which illustrate the 
general method in section 2. In section 3 we consider generalizations and further conceptual 
discussions of the general method. In section 4 we consider further examples. Section 4 can 
be read independently of section 3. 

The energy is only one out of the infinitely many charges that characterize the flux tube 
states. The other charges can be encoded in the generating functions of conserved charges 
constructed from the so called monodromy matrices. At strong coupling these monodromies 
are the main ingredient of the Y-system solution [12]. More precisely, at strong coupling 
one computes the holonomies of the flat connection between edges of the polygon, see figure 
2. In order to connect the weak and strong coupling descriptions, we study these objects at 
weak coupling in sections 5 and 6. In section 5 we identify these holonomies between edges 
at weak coupling and compute them for the simplest Wilson loops. The result we find is 
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strikingly different from the strong coupling analogous holonomies! This is not surprising. 
As mentioned above, at weak coupling we have a small finite number of excitations flowing 
from the bottom to the top whereas at strong coupling we have densities describing the 
propagation of infinitely many particles. To bring the weak and strong coupling descriptions 
closer - and hence provide further support for our identifications of section 5 - we consider 
in section 6 a scaling limit where infinitely many particles flow already at weak coupling. In 
this limit we insert, by hand, a large number of insertions in the loop. This would be the 
kind of situation one would encounter when studying N'^MHV amplitudes with /c ^ 1 and 
is the analogue of the Frolov-Tseytlin limit [13] in the spectrum problem. In section 6 we 
compute the weak and strong coupling holonomies in this limit and find a match between 
the two. 

In a very nice paper [14], Belitsky explored the connection of Null Polygon Wilson loops to 
the correlation function of two Wilson lines from the OPE point of view (we will also use this). 
He focused on MHV amplitudes with single particle excitations whereas we will consider 
mainly multi-particle excitations for N'^MHV amplitudes. He borrowed the technology of 
renormalization of Wilson lines while we will mostly explore the connection to Integrability 
for local operators. Both approaches are obviously tightly related and it is definitely worth 
exploring this connection further. 



2 From Wilson loops to two point functions 

We start with a correlation function between two Wilson loops in the fundamental represen- 
tation with two adjoint field insertions'^ 

Wbot = Tr 0{xi)W[xi,X2]0{x2)W[xi,X2]^ and Wtop = Tr 0\yi)W[yi,y2]0\y2)W[yi,y2]^ 

The Wilson line [xi, 0:2] connects points xi and X2 along some path. The conjugate Wilson 
line W^[a;i, ^2]^ connects X2 and xi along the same (but reversed) path. In the limit where the 
points Xi and yi become light-like separated, the correlation function (WbotVVtop) develops 
two poles corresponding to two fast particles going between C(xj) and O'^iiji). Each particle 
is charged under the gauge group and its interaction with the gauge field is approximated 
by a pair of Wilson lines in the fundamental and anti-fundamental between Xi and yi [15]. ^ 
This joins Vr[xi,X2] and Vr[?/i,?/2] into a single closed loop W and Vr[xi,X2]^ and Vr[yi,?/2]^ 
into another identical closed loop W'^ . In the planar limit (VTl^^) = iyV)"^. In other words, 
we start from a configuration where Xi — yi and X2 — 1/2 are space-like. Then, the residue 
of the double pole when (xj — ?/j)^ — )■ is computed by the expectation value of one closed 
Wilson loop as (see figure 3) 

hm (>VtopWbot) ^ 2_ 

(WtopWbot) tree 

Here, (WtopWbot)trce stands for the tree level correlator. We take the Wilson lines W^[a;j,|/j] 
to be composed of a sequence of null edges. Then the resulting closed Wilson loop is a 



■^Alternatively, one can add massless probe quarks and consider the operators q{x)W[x,y\q{y). 
^We do not distinguish between a Wilson line in the adjoint and a pair of lines in the fundamental/anti- 
fundamental representations. 
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Figure 3: A relation between a correlation function of two Wilson lines with insertions and 
a polygonal Wilson loop. On the left we have a correlation function of two Wilson lines with 
insertions Wbot and Wtop- We choose these to be composed of a sequence of null lines. At 
their tips we have scalar insertions indicated by blue dots in the figure. We take a limit 
where each one of the two scalars in Wbot becomes null separated from a conjugate scalar in 
Wtop- In that null limit the correlator develops a double pole due to two fast particles going 
along the left and right edges in the figure, indicated by the dashed lines. The residue of the 
double pole singularity is computed by two copies of a polygon Wilson loop. The polygon 
loop is composed of the two open curves in Wbot and Wtop, closed by the right and left edges 
into a closed loop. 

null polygon. The relation (1) is a simple generalization of the Correlation function/ Wilson 
loop correspondence of [15], see also [14]. 

The expectation value of null polygon Wilson loops needs to be UV regularized. As 
explained in [15], depending on the choice of regularization, there may be some issues in 
taking the light-like limit. These come about due to a possible recoil of the fast particles. In 
this note we will always consider ratios of polygons such that these issues cancel out. 

We will now introduce the OPE vacuum. For that, we take the point xi to be null sepa- 
rated from X2 and the same for yi and 1/2, see figure 4. Then, the four points {xi, X2, yi, 1/2} 
define a null square called the reference square. By a conformal transformation, we fix these 
points to 

{xi,X2, yi, 1/2} = {(0, 0, 0, 0), (0, 1, 0, 0), (1, 1, 0, 0), (1, 0, 0, 0)} (2) 

where the first two components are lightcone directions denoted by + and — respectively. 
The Wilson loop along the reference square, denoted as W^^'^^^'^'^\ represents the vacuum 
state of the flux tube [2]. 

The bosonic polygon Wilson loop is dual to the MHV scattering amplitudes [17, 18]. 
The other non-MHV amplitudes are captured by a supersymmetric generalization of this 
loop [19, 20]. Formula (1) can be easily generalized for the super loop case. The new 
ingredient is that we will also have insertions of various fields at the edges and cusps of the 
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Figure 4: The setup of Wilson correlators considered in this paper, (b) On the right we 
have a square representing the OPE vacuum. In that case, the bottom and top operators, 
Wtop"^'^°'' and Wbot^'^'^^'' , are composed of a Wilson line along a single null line connecting two 
scalars. (a) On the left we replace that single null line by a sequence of null lines connecting 
the two cusps in the bottom or top of the square. 



top and bottom Wilson lines as we now describe.'' The super loop has components which 
are the dual of the several scattering amplitude components. That is, is a polynomial in 
dual grassmannian variables rjf where z = 1, . . . , n is a particle label and A = 1, 2, 3, 4 is an 
R-charge index. We should always have R-charge singlets so that the total power of r/'s in 
each monomial is a multiple of four. Then it is convenient to factor out the purely bosonic 
loop which is obtained when setting to zero all grassmannian variables. That is, 

w/Wuuv = 1 + vlv-vlvf n^''''^ + vlv-vMvlvWovt 7^(^^■'=')(™"°^) + . . . (3) 

The 7^'s are the so called ratio functions [22]. The first one is the NMHV ratio function, 
the second is the N^MHV one etc. These are finite conformal invariant quantities.^ If we 
furthermore divide by their tree level expressions 

and similarly for N^MHV and so on (4) 



(ijki) 

tree 



then the resulting functions r only depend on the conformal cross ratios.^ These are the 



^Insertions on the two null lines of the two fast particle can, in principle, be taken into account by 
adjusting the operators O at the ends of Wtop/bottom [15,21]. We will not consider these components in this 
note. 

^Alternatively, one could consider the BDS stripped amplitude [10] which is the supersymmetric general- 
ization of the Reminder function [18]. This is also a finite conformal invariant quantity. 

^Without dividing by the tree level quantity we would need to deal with Helicity factors as in [6]. In 
particular, the conformal transformation (6) introduced below should also act on these factors and so on, 
see [6] for details. 
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finite objects that we will consider to illustrate the OPE method. More precisely, we will 
consider a very particular set of components (i.e. choice of ijkl . . . ) which have a particularly 
clean OPE interpretation. Namely we consider the components multiplying the monomial 



n 



The indices ia are associated to the bottom and the indices ja belong to the top part of the 
polygon. To leading order, a pair vlaVfa+i corresponds to inserting a complex scalar = Z 
at the bottom cusp Pj^ while the pair V^a^ja+i corresponds to inserting the conjugate complex 
scalar (f)^^ = Z ai the cusp Pj^ in the top. The grassmannian variables associated to the flux 
tube edges are rji and they are never used. To leading order (at tree level), the 

scalars just propagate from the bottom to the top cusps respecting planar ity, see figure 5. 

For obvious reasons we denote such components as Z — )■ Z components. Below, we 
will consider the three Z — )■ Z components depicted in figure 5 as our examples. The 
generalization to other components poses no conceptual obstacle. 

To render the discussion simpler while retaining all the physics, we will from now on 
restrict our polygons to lie in an M^'^ subspace. The generalization to M^'^ poses no conceptual 
obstacle. The null polygons in M^'^ have their edges along two null directions: and . 
We use (2) so that the points xi and X2 are separated in the x~ direction, X2 and ?/2 are 
separated in the x^ direction and so on. 

As mentioned above, by a conformal transformation we set the first cusp of the bottom 
of the polygon to be at position Pi = (0, 0) and the first cusp of the top to be located 
at Pm+i = (1;1)- These coincide with two of the cusps of a reference square, see figure 
4. Then we generate a family of polygons by acting on the bottom cusps (only!) with 
the time translation symmetry of the reference square. This generates a family of polygons 
parametrized by T = e~^^ where the x~^ coordinates of the bottom cusps are given by 

4{T) = +"2^ + • (6) 

1 — xj + I xj 

This is a conformal transformation that leaves the points and 1 invariant. We now have 
a family of Wilson loops W{T). We will now study the ratio functions r[T). It admits an 
OPE decomposition. At leading order 

r(T) = 1 + / [log(T)D(T) + b{T)\ + 0{g') (7) 

The function D{T) is the one loop OPE discontinuity and arises due to the one loop correction 
to the energy of the flux tube excitations that propagate. To measure the energy of the 
excitation we should work out what the flux tube Hamiltonian is. This operator will act 
on the bottom or top Wilson lines which we now introduce. Since we are dealing with a 
ratio function there are two Wilson loops relevant for the ratio (3): (a) the supersymmetric 
Wilson loop and (b) the bosonic one or MHV loop. Accordingly there should also be two 

^However, for the OPE purpose we do not consider components involving i^^ and ?7m+i. That is we only 
consider components that do not involve insertions along the two null edges of the flux tube. 
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Figure 5: There is a simple class of N'^MHV amplitudes which we use as a laboratory for 
exploring the OPE physics. We denoted such components as Z ^ Z components, see (5). At 
tree level they are given by k insertions of scalars Z at the bottom cusps of the polygon and 
another k insertions of the complex conjugate scalars Z at the top cusps. At tree level they 
are given by exactly k free scalars flux tube excitations going from the bottom to the top 
part of the polygon. At one loop these particles start interacting and feeling the flux. Such 
components are pure examples of multiparticles in the OPE. Figure (a): NMHV octagon, 
(b): N^MHV octagon, (c): N^MHV dodecagon. The last two correspond to two particle 
examples. These three examples will be discussed in detail in this paper and illustrate the 
general method. 



types of Wilson lines in both the top and the bottom. To leading order in the coupling these 
Wilson loops read 



W, 



bot 



^MHV 



Tr [Z(0, 0) ^ Z«(T), a;r) . . . Z«(T), a;rj ^ Z(0, 1) ^] 
Tr[Z(0,0)*Z(0,l)*] 



(9) 



and similarly for the lines in the top. The * in these expressions indicates that the fields are 
connected by the fundamental Wilson lines along the bottom and top contours as depicted 
in the figures. 

At one loop, the Hamiltonian EI will act on each pair of neighboring Z insertions in these 
Wilson lines as EI = (7^ X]a=i ^a,a+i where L is the number of scalars Z on the bottom (or 
top) Wilson line. Hence we can focus on a single pair of neighboring insertions inside the 
trace. These can be either null separated or not. Lets first consider the case where they are 
null separated in the x'^ direction and omit the direction x~ that plays no role,^ 



Tr [. . . Z(0) * Z[x) . . .] = Tr 



E 1^(0)^+^(0) • • • 



0{9) 



(10) 



^Of course, the same discussion holds true with and x interchanged. 
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Figure 6: To compute any one loop OPE discontinuity one computes the expectation value 
of the Hamiltonian (H) = J2i=i{'^i,i+i) ■ For N^MHV amplitudes of the type described in 
the main text the length is L = k + 2. In this figure the component 7'(2>3,io,ii)(4,5,8,9) ^j^g 
N^MHV amplitudes for 12 gluons is represented. 



where is in the fundamental representation. These operators are very well studied. They 
can be mapped to the so called SL(2) subsector of the full A/" = 4 PSU(2,2|4) spin chain. 
The SL(2) chain is a non-compact chain where at each site we can have an arbitrary integer 
which is identified with the number of covariant derivatives. The corresponding Hamiltonian 
is discussed in appendix C.^° Once we act with "H on these two neighboring sites we can 
re-assemble the result back into a Wilson line easily. This is also reviewed in this appendix. 
At the end of the day we end up with 

1 

/dt 
— [Z(0) * Z{x) - tZ{0) * Z{tx) - (1 - t)Z{tx) Z{x)] (11) 
t(l t) 



This kernel coincide with the renormalization of a null Wilson line with insertions and is a 
well known result [14,24]. 

If the two insertions are not null separated then instead of (11) we will have now both 
kind of derivatives, and D_. The action of the Hamiltonian is much less studied but it 
can be read from the Harmonic action of Beisert [23] as discussed in the appendix C. In this 
case, we find 

no Z (0,0) i<Z{x,y) -- 

Z{0)^Z{x,y) 
+ [•••] 

^•^The energy we are interested in is the Twist which differs from the Dimensions by an integer, the Spin. 
Hence, since we are interested in the anomalous part, we can use the Dilatation operator of A/" = 4 SYM [23] 
to derive the flux tube Hamltonian. 



1 

'dt f dz r dw 



t T 2TTiz J 2mw 



(12) 



- Z{xt[l + zlyt[l + w]) ^ Z{x{l -t)[l + \y{l - t)[l + -]) 
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Here, the -k stands for a Wilson line connecting the two points. As we explain in the next 
section, the one loop result (12) does not depend on the path and one can drop the *'s from 
the right hand side of (12). The dots in the last line, [...], involve insertions of the Z's out 
of the M^'^ plane. It can be read easily from the Harmonic action [23] and do not contribute 
when we Wick contract the result with two Z's in M^'^, for more details see appendix C. As 
a result, it will not contribute in the examples we will consider. 

After acting with the Hamiltonian as in (12) we can Wick contract the result with two 
Z fields. We get 



^Zd^Z,^---^{noZa^Zi,)^...) (x+ - 1) log X+ - (X" - 1) log X" I Q^^2) (^3) 



-kZd-kZ^-k-'-'kZa-kZb'k...) 

II II' 



where Za = Z{x^, ) and so on. Finally is the cross ratio 

{^t - ^t){^t - ^ir ^ ^ 

and similar for ■ The result (13) will be our fundamental building block. The Hamiltonian 
is Hermitian so we could have acted on the pair Zd -k Zc instead. If the two points in the 
bottom (or in the top) are null separated the result (13) simplifies considerably. For example, 
if x'^ = x'l we have = 1 and (13) reduces to the single logarithm logx~- This last result 
can of course be derived directly from the simpler representation (11). 

The OPE discontinuity in (7) comes from the correction to the energy of the bottom 
(or top) operators. We are computing a ratio function which means we should compute the 
average of EI = Yl!i=i in the super loop and subtract the same quantity in the bosonic 

loop. This leads to 

(top|e|bot) (top|e|bot)MHV 

Dil ) = — — ^ — ^ . (15) 

(top I Dot) (top|bot)MHV 

To compute the first average we sum over the action of "H on all nearest neighbor scalars in 
the bottom of the super Wilson loop as depicted in figure 6. The second average is computed 
similarly using the bosonic loop. For the action on each such pair we use (13). 

It is now very simple to compute the OPE discontinuity of any component with an 
arbitrary number of Z fields in the bottom. The Hamiltonian density EI acts locally on a 
pair of insertions at a time and the total number of insertions is irrelevant. This means we 
now have control over multiparticle states. The first interesting multiparticle states appear 
for the N^MHV amplitudes at tree level. We can now easily promote them to loop level. 

For example, for the N^MHV component of the 12 gluon amplitude de- 

picted in figure 7, we get the one loop discontinuity 

D{T) = = — f + log (Xl346X236lXl543j ) (16) 

X2365 X1254 

This is worked out more carefully in section 4.3. The result (16) agrees precisely with the 
unpublished results by Simon Caron-Huot and Song He! 
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Figure 7: The Dodecagon Wilson Loop. We generate a family of polygons with twelve edges 
W{T) by acting with a conformal transformation (twist) on the bottom part of the polygon 
(6). That conformal transformation leaves the square invariant and correspond to the "time" 
translation of the OPE. By a global conformal transformation we always set the reference 
square to be given by the cusps (0, 0),(0, 1),(1, 1) and (1,0). In particular this fixes the 
position of two of the cusps of the polygon, in this figure cusps 12 and 78. In the main 
text we considered the component 7-(2.3,io,ii)(4,5,8,9)_ j^^ leading order (tree level) this N^MHV 
Wilson loop is given by the propagation of the two scalars from the two cusps in the bottom 
to those two in the top. These scalars are denoted by the blue dots in this figure. 

Following [6], it should be simple to fully bootstrap all N^MHV amplitudes at one loop 
by considering the OPE expansion in several channels and making use of the SUSY Ward 
identities [25]. We can also work out the log(T)^ discontinuity at two loop level (and so 
on at higher loops) by acting twice (and more times) with the one loop Hamiltonian on 
the bottom of the loop. Later we will consider one such two loop example for an octagon 
N^MHV component. These will be new predictions. 



3 Generalizations and Comments 



In the previous section we have used the known result for the PSU(2,2|4) spin chain Hamil- 
tonian acting on the single trace operators to read the flux tube Hamiltonian acting on the 
bottom or top Wilson loops. In this section we will elaborate on the details of the map 
between the Wilson loop and the single trace operators. We will start from the vacuum 

(square) _ ^-^^^ 

move to the full polygon loop W, expressed as a sum of excita- 
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tions on top of the vacuum. The reader who is most interested in further apphcations and 
examples can skip this section and move directly to section 4 in a first reading. 

The reference square representing the flux tube vacuum is parametrized by the four points 
(2). As in (1) we have 

/ VW ('^'5'^'^''*') Vy (square) v 

lini , , ' = ))2 (17) 

The square bottom and top operators are simple null Wilson lines given by 

^(s^quare) _ [Ze^^-'^""- Z] = ^ ^ / Tr[ZD^Z] 



ml 

m=0 



oo 



(0,0) 



^(^quare) _ [Ze'^^-^)^- Z] = V ^^J^ Tr 



m=0 



(1,1) 



where we have Taylor expanded the Wilson loop operators in terms of single trace operators 
inserted at Xi = (0,0) and 1/2 = (1,1)- That is, in (18), we have mapped the expectation 
value of the Wilson loop W^''^''''"'^ into a sum of two point functions of local operators. 

Here, e = {xi - ytf is the null regulator. That is, (Wt^^;^"^''^>v£f ''''^)tree = ^ + 0(1). In 
(18) we have an arbitrary number of D_ derivatives acting on the second operator. However, 
the pole in the correlator (1) is dominated by the terms with infinite number of derivatives, 
displacing O from xi to X2- These are exactly the large spin operators with finite twist that 
are dual to the flux tube.^^ 

Suppose we now add one scalar Z insertion at a;~ = x on the bottom of the square and 
one Z insertion aX x~ = y on the top as prescribe in figure 8. That is, we now consider the 
bottom and top operators 



m! k\ 

m=0 k=0 



(19) 



(0,0) 



m=0 k=0 



(1,1) 



We now have 



/Txr\ 1- (VVtopWbot/tree 1 /^nN 

(W^)tree = T ^ ^ = (20) 

^ ^_>0 /-|/\;(square),,, (square) V V — X 

X'^^^top '^^bot /tree » 

That scalar propagator (20) represents the free propagation of the scalar Z excitation from 
the bottom to the top. The flux tube energy of that excitation is the twist of the corre- 
sponding operator Wbot or Wtop, minus that of the vacuum Wsquare- At tree level, it is just 



^^The map (18) also give us an alternative definition of the null correlator (17) (and (1)). Instead of taking 
the null limit e = {xi — — > 0, we can regulate the null correlators by putting a cutoff on the spin S of 
the operators in these sums. Then, the double pole divergence ^ is replaced by an S'^ divergence. When 
constructing UV safe ratios all these divergences cancel out leaving behind regulator independent quantities. 
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Figure 8: The square loop obtained as a correlation function of W^^p^^^^^^ and yvj^^^^^°\ We 
add a scalar ^(x) on the bottom line and a conjugate scalar Z{y) on the top. At tree level, 
these two scalar insertions are connected by a free propagator l/{y — x). The null correlator 
(^top^'^'^'^^^bot"'^'^'^^) is regulated by placing the bottom right scalar at a;" = 1 — e and the top 
left scalar at a;~ = e. 



equal toE = A — S=l + 0{g^). The spin cannot get any loop corrections. Therefore, 
loop corrections to the energy of the flux tube excitations are the same as the anomalous 
dimension of the corresponding operator. 

To compute the one loop correction to (W) due to the energy of the flux tube state, it 
is convenient to represent the operators Wbot and Wtop as spin chain states. We will denote 
these by (top| and |bot). The standard representation for the single trace operators in (19) 
as SL(2) spin chain states simply counts the number of derivatives 



,nL) = 



rTr 



ni! 



nil 



D"'Z...D1'^Z 



(21) 



So we have, for example 



|bot) = J2 (l-e)"'x'=|0,A;,m) 



(top|= (e-ir(l-#K^,0| (22) 



m,k=0 



m,k=0 



Then, the Wick contraction of the fields between two operators is represented by the overlap 
of the corresponding spin chain states, for more details see appendix C.^^ The correction to 



^^The states in (21) would have been an orthonornial basis provided that the two operators were inserted 
at zero and infinity. However, we have chosen to insert the two operators at zero and one. The two choices 
are related by a conformal transformation that relates the two frames. The corresponding orthonormal base 
in our (1,0) frame reads 

n 
k 



|n)) = e^+Hn) = E 



fc=0 



|A:) , ((m|n)) = 5„ 



(23) 



The Hamiltonian commute with the SL(2) generator L+i and is blind to such a change of basis. 
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(W) due to the one loop energy of the states is given by (top|EI o bot) where EI is the one 
loop spin chain Hamiltonian. In that way, we obtain the kernels (11) and (12). 

The one loop Hamiltonian kernel (12) acts on a pair of insertions that are separated both 
in the x+ and x~ directions. As stated before, it does not depend on the path connecting 
the insertion (parametrized by the 'V in (12)). That sounds strange at first since the 
gauge field in the Wilson line connecting the two points does contribute to the one loop 
Hamiltonian. To better understand what is going on, consider one such pair of insertions 
connected by a Wilson line Z{0, 0)W[{0, 0), (x, y)]Z{x, y). We first Taylor expand the Wilson 
line in covariant derivatives acting at (0, 0). Different paths correspond to different ordering 
of the covariant derivatives at (0, 0). For example, one of the operators in the expansion can 
be 

ZD+D^Z = ]^Z{D+,D^}Z + \gZF+^Z (24) 

Now the main point is that the two operators on the right hand side cannot be mixed by the 
one loop Hamiltonian! Even though they carry the same charges, they transform in different 
representations. The first is represented by two sites of a spin chain while the second is 
represented by three sites. At one loop, the spin chain Hamiltonian does not change the 
length of the chain (this does happen at higher loops [26]). As such, the second term with 
the F+_ insertion can only contribute at higher loops and can therefore be dropped. The 
first term, with the derivatives symmetrized, is independent of the order of derivatives we 
started with on the left hand side. The same argument applies for any number of derivatives 
and hence the result is independent of the path. 

In the previous section, we concentrated on computing the leading OPE discontinuity for 
non-MHV amplitudes. These are simple as they do not depend on the shape of the loop as 
explained above. The shape of the loop is important when considering the sub leading OPE 
discontinuities at higher loops or already at leading order when considering MHV amplitudes. 
For describing these we will now consider the bosonic loop dual to MHV amplitudes. 

Let us concentrate on the operator at the bottom, Wbot- Suppose that instead of the single 
edge on the bottom of the square we consider a bottom with four edges {^i, ^2, ^3, ^4} = 
{(a, 0), (0, 6), (—a, 0), (0, 1 — 6)}, see figure 9. a. The corresponding bottom operator is located 
at (0, 0) and reads 

Wbot - W^r"^^^ =Tr [Ze-^^+e^^-e-^^+e^i-^^^-Z] - Tr [Ze^-Z] (25) 

+ ... 



00 h. 
a" 



m=0 n=l k=l 



where we have subtracted from Wbot the vacuum yvj^^^^^°\ In the second line we have 
expanded the difference to leading order in the number of excitations. I.e., to leading order 
in the number of the corresponding spin chain sites. That leading term in (25) is basically 

"Abelian". It is equal to the integral of on the square of sides a and b as dictated 

by the Stokes theorem. Once contracted with the top, this term correspond to a single 
particle insertion on top of the flux. Note that the covariant derivatives in D^^^ D'^^ F^_ 
are symmetrized. They correspond to a descendant of the primary F+_ represented on a 
single site of the spin chain. Finally, the two terms in the commutator in (25) correspond to 
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Figure 9: Bottom Wilson lines. The two endpoints of the bottom and top Wilson lines 
are null separated from two other points in the top part of the polygon. The dashed lines 
represent two propagators that connect the two bottom endpoints to the two top endpoints. 
The two bottom endpoints (and similar for the top) are also null separated. However we 
choose them to be charged in such a way that there is no propagation from the left to the 
right. The simplest way to do this is to set the bottom insertions to be the complex Z scalars 
and the top insertions to be their complex conjugate Z. Note that the bottom part of the 
Wilson loop is made of the edges of the bottom plus a finite part of one of the two null edges 
that source the flux tube. That is, k\ in figure (b) is not the full edge. 



the insertions on the front and back copies of the bottom Wilson line, as indicated by the 
double line in the figure 9. 

At the next order we have operators with two excitations 

Wbot Itwo Excitations = / E ^"/^^^ ^r Z {T^^ [^^^ 6^" Z\ ] , (26) 

cj/37(5 

where J-""^ = b^'a^ D'^^^ D'^^ F^^ . The coefficents Cais-yS are obtained from (25) by symmetriz- 
ing the derivatives in the (+) and (— ) directions, keeping only the terms with two F+_'s. 
The general expression can be neatly expressed as a sum over paths with integration over 
a surface bounded by the path, but it does not seem very illuminating.^'^ Once contracted 
with the top, the operator (26) correspond to two flux tube excitations. 

Let us finish this section by enumerating the three types of loop corrections in our ap- 
proach entering the Ratio function (4) or its bosonic counterpart, the Reminder function. 

1. Corrections to the energies of the flux tube excitations. These are obtained by 
acting with the flux tube Hamiltonian EI on the bottom or top states. At one loop, EI 
only acts on the nearest neighbors but at higher loops the range of interaction grows. 

2. Bare form factor. These are obtained by first Taylor expanding the bottom and top 
Wilson operators. For example, for the bosonic bottom Wilson operator parametrized 

For the reference we quote only the first few terms in the expansion (26): cim = 3/8, C1121 = 1/4, 
C2111 = 1/9, C1112 = 1/6, C1211 = 2/9 etc. 
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by an ordered set of edges^'' ki,k2, ■ ■ ■ ,ks (see figure 9.b), that expansion takes tlie 
form 



\, {ki-Dr^ {ks ■ DT' ^ 



mil rris 



(27) 



Wbot = Tr [Ze^'-^e^^-^ . . . e^'^^'^Z] = J2 " ' J2 

mi=0 ms=0 

Then, by symmetrizing the covariant derivatives we expand (27) in the number of flux 
tube excitations or spin chain sites. This expansion includes contributions at any loop 
order. 

Renormalization of the bare form factors. The local single trace operators ob- 
tained in the Taylor expansion of the Wilson operators (25), (27) are not normal 
ordered. As a result, the decomposition of the bottom and top in flux tube states also 
gets quantum corrections due to self contractions of the fields in the trace. We hope 
to come back to this type of corrections in the future. 



4 Examples in Greater Detail 

At one loop order, the ratio function can be decomposed as in (7). Similarly, at any loop 
order, we have 

oo 21 ' 
1=0 ' m=0 

where D„i(T) admits a regular expansion around T = 0. This follows from the general 
arguments of [1] generalized to the supersymmetric loops, see e.g. section 5 in [6]. At each 
loop order the m = / term is called the maximal OPE discontinuity and can be computed 
easily. It is simply given by 

(top|e-^^^-«|bot) (toplbot) ^ ^ [-2Tg^]' / x 

(t0p|e-2.^-^|b0t)MHV^t0p|b0t)MHV l\ I ^ + ""^d } ) [^^ ) 

For example, 

^(1) ^ (top|e|bot) (top|e|bot)MHV 

^ (top|bot) (top|bot)MHV 

^(2) ^ (top|e^|bot) (top|e^|bot)MHV ^^(1) (top|e|bot)MHV .3^. 

^ (top|bot) (top|bot)MHV ^ (top|bot)MHV 

These averages can all be easily computed as explained after (15) and elaborated further 
below. As already emphasized above, a great advantage of writing things as in (29) is 
that the single particle and the multiparticle states are treated on the same footing. We 
have already worked out one example of a one loop discontinuity d\^^ above, see (16). In 
this section we will provide more details on that computation and also present a few other 
examples which will illustrate a few novel structures. 



^^A small comment: Since we have fixed xi and X2 to be light-like separated, the first (or the last) edge 
of the bottom is not the full edge of the original polygon but only part of it, see figure 9; similarly for the 
top operator. 
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Figure 10: Left: The octagon in M}'^ kinematics. The shaded region is the reference square. 
In the a;"*" direction we used the SL{2)^ subgroup of the conformal group to fix three of the 
points. The last one is parametrized as in (6). Right: The r^^^^^^ component at tree level. 
At tree level is corresponds to a scalar insertion at the cusp 23 that propagates from the 
bottom to the top cusp 67. There are also two other scalars at the endpoints of the bottom 
and top Wilson lines that generate the left and right fast particles that complete the two 
Wilson lines into a closed Wilson loop; these fast particles source the flux tube on top of 
which the middle scalar propagates. 



4.1 NMHV Octagon - Single Particle. 

As a first example we consider the octagon NMHV ratio function r*^^^^^^. For NMHV ampli- 
tudes we have, at leading order, a single scalar propagating from bottom to top, see figure 
10b. In figure 10a we have used the 5*^(2)+ conformal symmetry to fix three of the x"*" 
coordinates. The remaining cross-ratio is parametrized by the fourth coordinate which is 
given in terms of T as discussed around (6). For this example, we can choose the reference 
square as depicted in figure 10. 



4.1.1 One loop 

All the insertions are located at the bottom and top of the reference square. As explained 
before, in this case, (13) reduces to a single logarithm 



( ...-k Zd-k Z^-k ■ ■ ■ -k {H o Zg-k Zi,) -k ... ) _ {Xg - - x^ ) 

{...-k Zd-k Z^-k ■ ■ ■ -k Za-k Zb-k ...) {x- -x-){x~^ -x^) ' 
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The two insertions at the endpoints of the bottom Wilson hne are at x~ = and x~ = 1. 
Similarly, the two insertions at the top are also at a;~ = and x~ = 1; they are null separated 
from the bottom insertions. Hence, when acting with the Hamiltonian as in (32) we will get 
UV divergent results. These are physical divergences due to a divergent anomalous energy of 
the flux tube vacuum. When we construct the ratio function we subtract off that energy and 
should end up with a finite conformal invariant quantity. To see explicitly how this comes 
about we regulate these divergences by shifting slightly the position of the top and bottom 
insertions. We put the two insertions at the endpoints of the bottom Wilson line at x~ = 
and x~ = 1 — e and the insertions in the top at x~ = 1 and x" = e, see figure 8. 

From (30) and (32), one can immediately read off the two-point functions of the top and 
bottom Wilson lines in the presence of the Hamiltonian insertions, 

(H) = (Hi,2) + (7/2,3) + (Hs.l) 

xj(a£-e) (1 -x^)(l -xj -e) ^ 2e - 1 
= log . - _ +log f - _ -\ +^°g 2~' (33) 

vacuum = (EI)mHV = ('^1,3)mHV + ('^3,i)mHV = 2('Hi_3)mHV = 2 log — — (34) 



where 



_ (top|e|bot)MHV \ _ (top|'Hi,2|bot) 

= (toplbot)MHV ' ^^'''^ = (toplbot) ' 

and so on. As anticipated, we see that all epsilon divergences neatly cancel out when we 
subtract (34) from (33). At the end of the day, the one loop OPE discontinuity (30) for this 
component is given by 

= log [t^^I • (36) 

That is the coefficient of logT = logXi432 in the one loop result. This agrees perfectly with 
the well known results for the NMHV ratio functions, see e.g. [27]. 



4.1.2 Two Loops 

The / loop maximal OPE discontinuity is obtained by simply acting with the Hamiltonian / 
times, see (29). For example, at two loops we need to deal with (H^) = Xlj j('^^,^+l'^jJ+l)■ 
If i and j are very separated then the two Hamiltonian densities do not communicate with 
each other and we can simply use the one loop results twice. If i and j are close there are 
two different cases: (a) the two Hamiltonians are on top of each other (for i = j) or (b) the 
two are not aligned (for i = j ±1), see figure 11. 

For the NMHV component r^"^^^"^^ there is one further simplification: All bottom insertions 
are located at the same x~^, see figure 10. The same is true for the top insertions. In this 
case we can use (11) and (32) to readily obtain the new building blocks mentioned in the 
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(a) (b) 

Figure 11: At two loops one needs to compute the expectation value of a product of two 
Hamiltonian insertions. The two nontrivial cases are depicted in this figure. In figure (a) we 
depict the case {T-Qij^i) where the two insertions sit on top of each other. In figure (b) we 
have the case (7{j^j+i'Hi+i,i+2) + (^i+i,i+2^i,i+i) where they are slightly misaligned. 



previous paragraph, (see appendix B for more details) 



(^... -k Z^-k Zc-k ■ ■ ■ -k Za-k Zi)-k ...'^ 



-2Li2(l-x 



abed J ' 



(37) 



and likewise, 



r 



1 



_ I 

{ ...Zf-k Ze-k Zd- . . {T-Lb,c O T-igfi O Zg-kZb-k Zc)... ) 

(...Zf'kZe*Zd...Za*Zb'k Zc-..) 

^ I I I I I I ^ 



log Xabde log Xbeef + ^'^^{l " Xabef) (38) 



where the cross ratios are defined as in (14). Now, to compute the OPE discontinuities 
at higher loops we need to sum over the average of all the possible combinations of the 
Hamiltonian insertions. At two loops, the following terms involved in (31) are 



j+2,j+3j 



(h2)mhv = mh) 



MHV 



(39) 



Putting together all the ingredients, one obtains the following expression for the maximal 
two loop OPE discontinuity (31) 



D 



(2) 



log' 



^ Xl243 

(x 



1243 J 



1 2 

-log'(l-X^243) + 371■^ 



(40) 



In appendix A we cross check (36) and (40) against the usual OPE promotion. It would be 
interesting to cross-check this prediction against (more) direct computations. 



4.2 N^MHV Octagon. Two particles. 

We now move on to the 2-particle states. We consider the component (^72^3^7%)('73'74'76'77) 
of the N^MHV 8-point amplitudes, see figure 5. 
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4.2.1 One Loop 



The logic of section 4.1.1 goes through without any modification for this case. The fact that 
we are now dealing with two particles instead of one makes no difference since all we have to 
do is consider the local insertions of the Hamiltonian densities; they do not care about the 
total number of insertions. The 1-loop OPE discontinuity obtained in this way reads 

(Xl243 ~ l)Xl243 



D^^ = log 



(X 



1243) 



(41) 



4.2.2 Two Loops 

We now move to two loops. At two loops we also need to compute ('Hj^j+i'Hjj+i). Two such 
examples are depicted in figure 12. For the octagon example under consideration, it turns out 
that in all cases that one needs to consider, one of the Hamiltonian insertions always acts on 
two insertions which are null separated. For example, in figure 12a the bottom Hamiltonian 
acts on the two insertions located at the same while in figure 12b it acts on the two at the 
same x~ . Hence to compute (T/j^j+iT^jj+i) we replace the action of one of the Hamiltonian 
insertions on the two null separated points by (11). This leads to an integration over t. Then 
we use (13) to get rid of the second Hamiltonian and finally we do the integral over For 
example. 



(top|?{l,2?^2,3|bot) 

(top|bot) 



Lio 



1 - 



Xa I Xn 



I Xr) Xa 



- Li2 



l Xi243 "; 



Xa I Xn 



[Xo Xa }€- 



1 



log'x^243 (42) 



where e is a regulator, see section 4.1.1. When we add up all contributions and compute (31) 
it drops out. All other contributions can be easily computed, see Appendix B for some useful 
building blocks. Summing up all the contributions as in (31), the 2-loop OPE discontinuity 
for this component is predicted to take the following form. 



D. 



(2) 



log^(l - X1243) - 2 log 



1 X1243 



logXi243 + logXm3 log 



X1243 



Li2(l-X^243)-(43) 



2 X1243 (Xl243)^ 

In principle, there is no conceptual obstacle to generate in this way infinitely many higher 
loop predictions. 



4.3 N^MHV Dodecagon 

Let us close this section by giving one final example, the component 7-(2:3,io,ii){4,5,8,9) ^^le 
N^MHV amplitude for 12 particles at one loop, see figure 7. The computation is basically 
the same as the previous two examples. But since now the scalar insertions are separated 
diagonally in both the top and bottom, the rational functions multiplying the logarithms in 
(13) survive and the OPE discontinuity at one loop leads to (16). 

^^Note that if neither the top nor the bottom insertions were null separated this computation would be 
a bit more complicated. More precisely, after acting with the first Hamiltonian we would get two insertions 
which might leave the M.^'-^ plane. Then, the dots in the last line of (12) will also contribute. For all examples 
we will consider we can avoid dealing with them as explained in the text, see also section C.2.1. It would 
certainly be interesting to consider an example where these terms need to be taken into account. 
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Figure 12: Examples of Hamiltonian insertions for the octagon N^MHV Wilson Loop at two 
loops. To compute the two loop OPE discontinuity we should sum over all possible pairs of 
insertions. For the octagon all cases exhibit an important simplification: at least one of the 
two Hamiltonian insertions acts on a pair of null separated points. For example, in figure 
(a) the bottom insertion acts on two points null separated along x~ while in figure (b) the 
bottom insertion acts on two points null separated along x~^. Because of this nice feature we 
can always use (11) followed by (13) and easily compute all building blocks, see appendix B 
for more details. 

5 The Monodromy Matrix 

We have seen that it can be quite useful to think of the Wilson loop expectation value as a 
correlation function of two Wilson lines, a bottom and a top one. These Wilson lines can be 
expanded in local operators. These local operators can be thought of as non-compact spin 
chains and a great deal is known about these from Integrability [28]. For example, above 
we used the knowledge of the action of the Dilatation operator on these spin chains to learn 
about higher loop predictions for Wilson loops. In particular we saw how to tame multi- 
particle corrections in the OPE. The Dilation operator acts as a spin chain Hamiltonian on 
the local operators and measures the energy of the state produced at the bottom and flowing 
to the top. This is one of the conserved charges of the spin chain. The = 4 spin chains 
are integrable and there are many other interesting charges that are well understood and 
that one might envisage measuring and making use of. In Integrable models, the generating 
functions of such higher charges are the so called monodromy matrices and transfer matrices. 
In this section we shall initiate the study of these objects for Null Polygonal Wilson loop at 
weak coupling (to leading order in the 't Hooft coupling). 

For guidance and further motivation, it is very instructive to recall how the story goes at 
strong coupling where the full artillery of Integrability was put to use [5, 12,29]. For large 
't Hooft coupling, the Wilson loop is governed by a minimal surface area in AdS [17]. The 
main ingredients entering the solution to the strong coupling problem are holonomies of the 
fiat connection on the worldsheet [12]. The world sheet has a disk topology and therefore 
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Figure 13: The bottom Wilson line along which k scalars Z are inserted at 



the holonomies are not closed. Instead, they can end on an edge of the polygon, see figure 2. 
When restricting to M^'^' kinematics, these holonomies go between two edges of the polygons 
in the same null direction. At each of the edges, they are contracted with a fiat section 
associated to that edge. That section is called "small solution" and is the unique fiat section 
that is decreasing as we approach that edge. To summarize, schematically, the fundamental 
building blocks are then the holonomies connecting edges i and j 



sl ■ £ 



Vexp A{u) -sr, SL = Si{oOi), sr 



(44) 



Here A{u) is the fiat connection used to transport the sections from one edge to the other. 
The connection depends on an arbitrary complex number u called the spectral parameter. 
Si is the small section associated with edge i; oOi is a point at edge i and e is the SL{2) 
metric Eajs = P — a with a, /3 = 1,2. Finally, the minimal area is some functional of these 
holonomies, see [12] for more details. 

We will now try to identify the analogue of the holonomies and small solutions at leading 
order at weak coupling. The hope is that once the strong and weak coupling result are 
expressed in terms of the same ingredients, we may understand how to extend them to any 
value of the coupling. 

To present the main ideas we will focus on some particularly simple components of the 
super Wilson loop where the states created at the bottom (and absorbed at the top) admit a 
very transparent spin chain description. At the end we comment on how to generalize these 
considerations to any possible components and even present some results for the bosonic 
Wilson loop. 



For simplicity, we start by considering the Z ^ Z components introduced in 



Fur- 
thermore we will consider the case where all the insertions at the bottom are located at the 
same x^ = and the same for all the top insertions, inserted at a:"*" = 1, see figure 13. In 
this case, to leading order at weak coupling, the bottom and top Wilson lines are simply^^ 



W, 



hot 



Tr [Z(xo, 0) * Z{xi, 0) * ■ ■ ■ ^ Z{xk, 0) * Z{xk+i, 0) 



(45) 



^^To be more precise, all the scalars Z are dressed by extra pre- factors that carry helicity weights [19,20]. 
However in practice, we will be considering ratios constructed from these Wilson operators, for example (52), 
and these multiplicative factors drop out. 
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where Xq = and Xk+i = 1 — e with e being a regulator that will drop out in any physical 
quantity. Equivalent ly, 



w,,t = y . . . i^kr^o_jy^r^^ izd^^z . . . di^zdi^+^z] (o, o) 

^ nil Ukl Uk+il 

nj=0 

and similarly for Wtop- For N'^MHV components, we are dealing with spin chains of length 
L = k + 2. In the OPE language, they correspond to k particle states in the OPE flux tube. 
In the spin chain notation, the bottom operator simply reads 

oo 

|bot) = [(^i)"' • • • (^-t+i)"'^'] |0' ^1' • • • ' ^fc+i) ' (46) 

nj=0 

see appendix C for more details. Such spin chains are the so called SL{2) spin chains where 
at each site we have a spin transforming in the spin s = — 1/2 representation of SL{2). The 
SL{2) generators act as La = X]n=i where lI"-* acts on site n. Acting on kets we have 

L_i|n) = {n + l)\n + l) L.ioZ{x,y) = Z{x,y) 

Lo \n) = in + l/2)\n) or Lq o Z{x,y) = {xd, + 1/2) Z{x,y) (47) 

L+i\n) = n\n — 1) L^ioZ{x,y) = {x^d^ + x) Z{x,y) 

when acting directly on fields. Of course, it is straightforward to go between these two 
representations. The action on kets is the convenient one if we want to use (46) while the 
action on the operators can be used to act directly on each of the fields in (45). 

The monodromy matrix can be thought of as the scattering of an auxiliary ghost particle 
with the physical excitations. It is given by 

L = ni{u)... niiu) , TZjiu) = Rj{u - i/2) (48) 

where u is the so called spectral parameter and where the i?-matrix Rj acts on a tensor 
product of two spaces: The physical space at the chain site j and an auxiliary space Vq 
associated to the auxiliary particle. The length L = 2fc + 2, see (45). The most important 
property of the i?-matrix is that it obeys the triangular relation known as Yang-Baxter. 
One can think about the i?- matrix as a discretized/quantum version of the strong coupling 
holonomy. Yang-Baxter is the quantum analogue of the fiatness condition at strong coupling. 

The monodromy matrix acts on the physical Hilbert space tensored with the space asso- 
ciated to an extra auxiliary particle. That is, denoting the indices of this auxiliary space by 
capital letters A,B, . . . and the indices of the physical sites by ia, ja, we have 

(^n:/^)I= E (^lW)A■^(^2(«))^^..(7^.(n));^^,^ (49) 

Ai,...,Al 

The ghost particle transforms in some representation of SL(2). Different representations 
correspond to different auxiliary spaces Vq and define different monodromy matrices. Which 
one is the relevant one for our purposes? There are two obvious candidates that deserve to 
be analyzed. 
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One natural choice is to take the representation of the auxihary space to be the same 
as in the physical space. That is, the non-compact, unitary, spin —1/2 representation. 
In this case, A,B,... are of the same nature as the physical indices ia and ja and take 
values A, B = 0,1,2, ... . The monodromy matrix in this representation is certainly a very 
important object. For example, the spin chain Hamiltonian and the next local higher charges 
can be obtained by expanding the (trace of the logarithm of the) monodromy matrix at a 
particular value of the spectral parameter. The i?-matrix for this choice of auxiliary space 
is given in [23]. It would be very interesting to study this case in greater detail. In this 
paper we focus on the other natural choice, which is mostly motivated by the story at strong 
coupling. 

The second natural choice is to let the auxiliary space Vq be the simplest /smallest it can 
be. Namely, we take the auxiliary particle to transform in the spin 1/2 compact (and hence 
non-unitary) representation of SL{2). In this case Vq = and the auxiliary particle can 
either be spin up or down, that is A, B = 1,2. At strong coupling, for polygons in M}'^ the 
small solutions were 2-component spinors which were transported between edges using the 
flat connection in this same compact representation. This is the main motivation for this 
choice. The spin 1/2 SL(2) R-matrix is very simple and is given by [32] 



=n{u) 



^u^ + 1/4 



^ j(aux) ^ j(phy) + 1 ^(aux) ^ ^(phy) 

O / J Oj Ob 



(50) 



where L^a"^^^^ are the conformal generators acting on the physical space as in (47). L^a^^^ are 
the conformal generators acting on the auxiliary space in the compact spin 1/2 representation. 
That is, they are simple Pauli matrices 

.(aux)_/0 l\ .(aux)_/^ \ . (aux) _ / 0\ , . 

Using (47) it is easy to read off the "expectation value" of the R-matrix after acting with 
it on a single scalar Z{x) in the bottom spin chain and contracting the result with the 
corresponding scalar Z{y) in the top spin chain. This way we get rid of the physical space 
and what is left over is simply a 2 x 2 matrix that acts on the auxiliary space only,^^ 



which yields 



(top|i?(^)|bot) _ Zjy) (nju) o Z{x)) 

(toplbot) " Z{y)Z{x) ^ ^ 



^ _ i x+y ixy 



^^For the spectrum problem the connection between the monodromy and transfer matrices in compact 
representations at weak and strong couphng was beautifully studied in [30,31]. 

""^^We omitted the dependence on the x+ direction which plays no role and drops out. 
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We can now construct the analogue of the most important building block in (44), namely 
the path ordered exponential integral of the fiat connection. It corresponds to contracting 
the monodromy matrix (48) with the top and bottom operators. These operators are made 
out of the several insertions at points Xi and Ui. That is, we simply need to multiply the 
several 'Kxi,yi{u) operators as 



{u) ■ {u)-...- M,,,,, {u) ■ R,,^„y,^, (u) = n{u). (54) 



Each of these (reduced) R-matrices propagates an auxiliary vector from one site to the next. 
In total this object carries a vector from the rightmost site to the first one. It can be thought 
as the analogue of the path ordered exponential in (44) used to carry the small solution at 
edge j all the way to edge i. 

Next we move to the next ingredient in (44), the left and right small solutions sl and 
Sr. At strong coupling when we are close to an edge and propagate a spinor all the way 
to that edge there are two possibilities. The spinor can be tuned in such a way that, once 
propagated all the way to the boundary, it vanishes. If it is not tuned it will explode instead. 
The directions which lead to the vanishing results close to the left and right edges are sl and 
sr in (44). We can now repeat this analysis at wealc coupling. The wealc coupling analogue 
of the propagation from a point close to the edge all the way to the edge is given by the 
action of the left-most and right-most M-matrices in (54). So, for example, we define sl [sr) 
as the spinor that is annihilated by the left-most (right-most) M-matrix once we propagate 
it from the right (left). 



» ■ SL = 0{e) , SR-e- R.,^,,y,^M) = 0{e) (55) 



Of course, these spinors are defined up to an arbitrary normalization. Such spinors exist 
because these matrices are singular as Uq—xq = Xk+i—yk+i = e — )■ 0. (Also at strong coupling, 
once we are propagating all the way towards the null edge we approach the boundary of AdS 
and the warp factor leads to a divergent contribution.) 

Using (53), we obtain 



These are the small solutions when the edges are located at a;o = and Xk+i = 1- We could 
repeat the exactly same computation for edges at arbitrary positions x. We find in that case 

X 



lj+0{e), (57) 

This result is a local property of the edge and is independent of whether we define it as 
acting from the right or from the left. The results for different a;'s are related by conformal 
transformations that relate the different a;'s.^° To study different OPE channels amounts to 
choosing different pairs of null edges. It would then be convenient to work with these more 
general small solutions. 

^^At strong coupling the background is classical and hence we never have to discuss the contraction with 
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Figure 14: Weak coupling transfer matrix (59) represented by the solid red line. The dashed 
(blue) line are the physical space lines while the solid (red) line represents the auxiliary 
space. This is the analogue of the strong coupling holonomies represented in figure 2. 

We now have all the ingredients to write down the analogue of (44) at weak coupling. The 
first and last M-matrices in (54) are used to define the directions sl and sr. The remaining 
M-matrices in (54) are the analogue of Ve^ ■^^'^^ in (44). That is, we have 



Of course, this depends on the overall normalization of the small solutions. We can construct 
a normalization independent quantity by dividing by their scalar product sl ■ £ ■ sr. This 
scalar product can be thought of as the overlap of sl and sr in the vacuum. That is, we end 
up with 



which we represented in figure 14. This is the object that has a nice strong coupling analogue. 
It would be wonderful if we could now extend this object from weak and strong coupling to 
any intermediate coupling. Also, we would need to know how to extract the Wilson loop 
expectation value from such holonomies, like we do at strong coupling, see discussion at the 
end. 

To get some better intuition/familiarity with this object we now present a few more 
examples and comments. In particular, in section 6, we will explore an example with many 
i?-matrices where we will be able to make direct contact with strong coupling. 

the top and bottom operators, they are automaticaUy done. Only the auxihary space is visible. Hence the 
reduced R-matrix M(m), which only acts on the auxihary space, is the analogue of the discretization of the 
strong coupling path ordered exponential. 

^"The conformal transformation relating the points reads: x — >■ ffqr^- This corresponds to the conformal 

transformation relating the small solutions s — )• g o s where g = ^ ^ • This specific representation of g 

follows from the fact that the physical and auxiliary space are entangled as in (50). For a related discussion 
see section 5.2. 



(L, -R)wcak = Sl- £ ■ ^x^,yA'^) " ^X2,y2{'^) " " " ^Xk,yk{'^) " ^R 




nu) 



Sl-£ ■ Sr 



(59) 
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5.1 The Simplest Examples 



The simplest possible example is the vacuum. That is quite uninteresting. It corresponds 
to putting no insertions at all so that in (59) we have no M-matrices and hence T{u) = 1. 
The next to simplest example is the case with a single insertion. Let us take the component 
V2V3V7V8 of the octagon depicted in figure 10 as an example. There is one scalar insertion 
between edge 1 and edge 5. Therefore here the monodromy matrix contains only a single 
]R('u) corresponding to this scalar. Then we have, 



1 



'^^M = ^ - + ^^1234 , (60) 

^yu^ + 1/4 L ^ 



where X1234 is the cross ratio defined previously. 



5.2 Gauge Transformations and Globally Defined Small Solutions 

The small solutions of the linear problem are the key ingredients in the computation of Null 
Polygon Wilson loops at strong couphng. In the previous section we identified the analogous 
definition of small solutions at weak coupling. More precisely, so far we only needed (and 
only computed) the small solution sj close to the corresponding edge j, see (56). At strong 
coupling, the small solution is defined close to the corresponding edge but can then be 
propagated anywhere in the worldsheet. Similarly, at weak coupling, we can identify the 
dependence of the small solutions on the position in the spin chain, as we now explain. 

As mentioned above, the analogue of the propagation with the fiat connection is given 
by the reduced M-matrices (53). Hence we can define the small solution s^'' at the nth site 
by propagating it k — n times from its value at the rightmost end of the chain, defined in 
(56). That is 

Sr^ = ^x.,yM ■ K..„+i,,„+,(«) ■ . • •K.„,,(m) ■ sn (61) 
In particular the constant spinor sr is the small solution for n = A; + 1. Similarly, for the 
left small solution we define s^"* = sl ■ £ ■ '^xi,yi{u) ■ ]Rx2,y2('") ' • • • I^x„_i,y„^i(^) ■ ^ ■ Then the 
weak coupling analogue of the brackets (L, R) as defined in (58) can also be written as 

(L,i?)^eak = si")-£-4"^ (62) 

exactly as at strong coupling. This quantity is independent of n. Similarly, the strong 
coupling brackets are Wronskians of the linear problem and as such they are independent 
of the world-sheet coordinates. The story is not yet complete. The M-matrices act on an 
auxiliary space defined at some given site and map that to a neighboring site. We are free to 
change our base for any such auxiliary space. In other words, (58) is invariant under gauge 
transformations 

^xuvi 9i-i ■ ■ , sr gk+\SR , sl^ Sl- gl (63) 

2 2 2 

where g^j^i is an SL(2) transformation (such that £ ■ gj^i ■ £ = g^^i)- In particular we get 



2 



(n) (n) (n) (n) —1 /r'A\ 

SR^9n-ys'ji\ Sl s\> ■£-g^\ (64) 
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which leaves (62) invariant. This is the direct analogue of the local gauge transformations 
at strong coupling which do not affect physical quantities such as well defined holonomies of 
the fiat connection. 



5.3 Open vs Folded 

Above we constructed an holonomy obtained by starting with a small solution at edge i, 
carrying it all the way to edge j and contracting the resulting vector with the small solution 
at edge j. In this procedure the two end-points play an important role. This is expected since 
we are dealing with open strings. On the other hand, we also know that for most practical 
purposes, a closed folded string behaves as two decoupled open strings. Basically, each side 
of the fold has an infinite effective length. Its two sides decouple since excitations cannot 
propagate from one side to the other in finite time. Therefore, it should be possible to express 
our main objects (59) in a language closer to the spectrum problem, where the key objects 
are closed holonomies. That is, traces of the monodromy matrix. This is the analogue of the 
usual adjoint versus fundamental story which we encountered already in the introduction. 
In the planar limit, these are trivially related. Indeed, it is possible to define (59) purely in 
terms of closed string like quantities, without ever introducing any small solutions as 

Note that in this expression we included the first and last M-matrices which are absent in 
(59); they depend on e through yo — Xq = yk+i — Xk+i = e. It is straightforward to prove 
(65).^^ That relation between the closed and open monodromies also makes clear that T[u) 
generates conserved charges. 

In the trace in (65) there are only R-matrices between Rxo,yo and Rxk+i,yk+i tiut none 
between Rxk+i,yk+i ^^owo- This is because, so far, we have chosen to put excitations 
only between Z(0, 0) and Z(l, 0), leaving the other side of the bottom and top Wilson loops 
empty, see for example figures 7 and 8. More generally, one could choose to have excitations 
on both sides of the bottom and top loops. In such cases, due to the projective property of 
Mi = '^xo,yo and Mr = '^Xk+-i^,yk+i) the trace of the monodromy matrix factorize as 



front back 

Tr ' ^ 



Tr [llfront(^^)] Tr [Vth!,ck{u)] 



Tr Ml ■ Mr] Tr [M^ ■ Mr] Tr 



(66) 



A symmetric choice - which we will make in section 6 - is to have a folded string where the 
back of the loop is just the mirror of its front. In that case, the front and back monodromy 
matrices are related as fiback(M) = ^front(— m)~^- 
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Here is a pedestrian explanation: Using (56), we have ■ e ■ M ■ sji — — Afn — M12 ~ — Tr \P ■ M] where 
the projector V — i \ ) . On the other hand, as e — we have IRa;;^^^,^^^^ [u) ■ M^x^.y^ (u) "^rrr x V. 



Hence (65) and (59) are equivalent; they can both be written as Tr \P ■ '&xi,yi [u) • . . . • Kajfc.yfc (u)] /Tr [7^]. It 
is easy to write the same proof in a way that does not rely on a specific choice of frame using the basis 
{s, £ • s} at each edge. 
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5.4 Monodromy for general Wilson loops and MHV amplitudes 



Above we saw how to define the analogue of the strong couphng holonomies at weak couphng. 
We did it for a particularly simple set of examples, namely we considered Z — )■ Z components 
of the super loop and, furthermore, we assumed that all the bottom excitations were created 
at the same (and similarly for the top). In fact, it is not hard to lift these simplifying 
assumptions as we now explain. 

First, suppose we still consider the Z ^ Z components in M}'^ but with insertions at 
arbitrary positions in the direction. Then the relevant i?- matrices would be SL{2)_ x 
5'L(2)_|_ i?- matrices which can be written as a four dimensional block diagonal matrix with 
each SL{2) occupying one block. Then the claim is that (59) is simply untouched. Namely 
the result only depends on the x~ coordinates of the insertions which are the variables 
appearing in (59). The reason is that the small solutions - associated to edges along the x"*" 
direction and localized in x~ - only have non-zero components in the SL{2)_ block which is 
the only block that becomes degenerate as e — )• 0. Hence the SL{2)^ block is simply never 
probed. This is related to the AdS^ reduction discussed in section 4.9 of [12]. 

Next we could consider different components of the super loop in M^'^ where the excita- 
tions are not the scalars Z but some of the other fields such as the fermions, field strengths, 
etc. The SL{2)cr -R-matrix still takes the same form as in (50). The only difference is 
that the SL{2) generators act on fields in other representations. These representations are 
parametrized by the conformal spin s. When acting on a field inserted at a cusp, instead of 
(47), they take the form 

Lj^i\n)s = {n + 2s -l)\n-l)s L+i o = [x'^d^ + 2sx] Os 

Lo\n)s = {n + s)\n)s or Lq o Os = [xd.^ + s]Os (67) 

L_i|n), = (n + l)|n + l), L_^oOs = d^O, 

where 

\n), = -\D10,) (68) 



Finally, as R is an SL{2) matrix, its overall normalization is modified from 1/^/1/4 + 



u 



2 



to 1/a/s(1 — s) +v?. Before, for the scalars, we had s = 1/2 but for the fermions s can 
be either 1 or 1/2, for example. Hence our previous computations can be straightforwardly 
generalized to these cases as well. When a field is integrated on an edge, the integration 
measure may shift the conformal spin of the state accordingly. For example, for the one loop 
MHV octagon we have a -D^F+_ insertion integrated along the bottom of the square. That 
field has conformal spin s = 1, but the integration shifts it back to s = 0. In that case, using 
the expansion of the bottom state (25) and similarly for the top, we find 



I e ill 



M (1 - e ^'^) log(l - e ^'^) tmhy V 2u J 
where a = — ^logxr342 is the cross ratio in the x~ direction (see figure 10) and tmhy = 



— log(l + e ^'^) log(l — e ^'^) is defined as the following logarithm of a ratio of polygons. 
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(see [1,4] for more details) 



r-MHV = log( 



)-log( 



(70) 



It measures what is flowing from the bottom to the top for the case of the bosonic loop. 
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6 Many Particles and Weak/Strong Coupling Match 



At this point let us summarize the simplest non-trivial holonomies encountered thus far. At 
weak coupling the simplest holonomies are those that measure single particle states. These 
are the tree level NMHV octagon (60) and the one loop MHV octagon (69). These expressions 
are both very simple. Up to a trivial multiplicative factor, they are given by a simple linear 
polynomial in u whose coefficients depend on the cross-ratios of the octagon. Let us now 
compare these weak coupling holonomies with a strong coupling one. At strong coupling, 
the simplest holonomies are those of the octagon MHV Wilson loop which (in some simple 
gauge) read [12,29] 



strong 



U 



exp 



dO' + e^°'^'^'^^'^'°s^i432-«smh(e')iogxi342 j 



2Txi 



cosh[^('u) - 6'] 



(71) 



where^^ 



d{u) 



1 u + 2g 

- log 

4 u-2g 



(72) 



Comparing this result with the simple polynomials obtained at weak coupling (69), we con- 
clude that the weak and strong coupling results could hardly look more different! 

We should not be surprised however. The simplest weak coupling examples describe the 
propagation of a single particle whereas at strong coupling we have an infinite number of 
excitations fiowing. They ought to look very different. 

In order to connect these two seemingly very different limits, and hopefully learn about 
the finite coupling interpolation, we would like to consider a limit where, already at weak 
coupling, we have an infinite distribution of OPE particles. For that purpose, we will now 
consider a scaling limit with infinitely many scalar insertions in the square Wilson loop.^^ In 
this setup, the weak and strong coupling holonomies will become very similar (even identical 
in some limit). This limit, which allows us to establish a bridge between weak and strong 

^^Note that the combination {u — |9cr) — {u — p/2) is the form of the spin chain transfer matrix for an 
excitation on top of the GKP vacuum [4]. There, p is the spin chain momenta of the hole. 

^•^In the strong couphng classical limit u is large (of order g) so that the spectral parameter which appears 
in the flat connection is of order one. 

^^Alternatively, we can have many Z insertions by going to very high MHV degree instead of just inserting 
these excitations by hand on the reference square. For the purpose of understanding the underlying physics 
both are fine but the latter is simpler and hence that is the example we will follow. 
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coupling is the direct analogue of the algebraic curve classical limit [30,31,33,34] which was 
of great importance in the study of the spectrum problem. 

We will consider the weak coupling side of this many particle example in subsection 6.1 
and the strong coupling side in section 6.2. As we will consider many complex scalar Z 
insertions, the dual string will have a non-negligible motion in the sphere. We will describe 
a simple string solution corresponding to a null polygon in AdS together with point-like 
motion around an equator in the sphere. We will see that in a particular limit, which is the 
analogue of the Frolov-Tseytlin limit [13], the weak and strong coupling holonomies match 
perfectly. 



6.1 Weak Coupling 

We would like to study a simple case where there are infinitely many particles flowing already 
at weak coupling. That is, we will insert infinitely many scalars Z on the bottom of the square 
and corresponding Z's at the top. There is a natural translation symmetry that preserves the 
two null edges that source the flux tube. Under this translation symmetry, points transform 
as 

^^^{o) = - _ , (73) 

1 — -I- e^'^Xj 

This is the analogue of (6) in the x~ direction. Note that this conformal transformation 
indeed preserves x~ = and = 1 which are the location of the two null edges. We would 
like to insert the Z's in such a way that preserves as much symmetry as possible. Therefore 
we consider insertions which are homogeneously distributed with respect to this a translation 
symmetry. More precisely, we will insert scalars Z in the bottom of the square at positions 

{a;~}bottom = < ^ ^ ^2a„ \ ^^^^ ^" = ^os{S)j , n = -J/2, . . . , J/2 , (74) 

Here log(S') is the cut-off in the a direction. At the top we consider insertions of scalars Z 
at positions 

where a is a free parameter, parametrizing how much the propagators from the bottom to 
the top are tilted. In (74) all points are inserted inside the interval (0,1). To make the 
configuration spacelike we consider the top part of the polygon to go from to 1 "from the 
outside". That is, in (75) we start close to x~ = at large negative a, go all the way to 
= oo at cr = and come back to x~ = 1 at large positive a. Note that from a conformal 



^^Since we are considering mostly polygons in M^'^, described at strong coupling by minimal surfaces in 
AdS^, the most closely related algebraic curve reference would be [34]. 

^^It is very much related to the spin cut-off discussed in footnote f 1 but not exactly the same as will 
hopefully become clear latter. We will come back to this point in a latter footnote. In practice this is a 
detail that does not play an important role in our discussion. 
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point of view there is no difference between the (75) and (74). Indeed they are related by a 
simple analytic continuation, 2cr — t- 2cr + + a. 

We want to consider the classical limit where the number of insertions is very large and 
so is u. This scaling limit was studied intensively in the spectrum problem, see e.g. [33] 
and [35]. For large u the reduced M-matrices (53) exponentiates 



x,y 



.U 



exp{iA^ei.k{u;x,y)), ^weak(M; a;, y) = - ""x^y (76) 



u 



x~y 2 x-y 



That is, if we have a smooth sequence of points (like we do) then the product (54) becomes 
a simple path ordered exponential integral, exactly as in (44)! The only difference is that 
the strong coupling connection A is now replaced by the weak coupling one, ^weak- This 
reinforces strongly the identifications of section 5. In sum, for our case we have simply 



h log{5) 



2 

ij j da A^ci^kiu; a) 



n{u) - ^x_j^2,y-j/2 ■'^x_jj2+i,y-j/2+i ■ ■ ■ ■■'^xj/2,yj/2 - ^exp 

-|log(5) 

Here, j = J / (log S) is the density of insertions. For example, at a = we have 

1 / -e2- 



(77) 



Xeak(«; ^) = ^ 2 sinh(2a) -e^^ ) " = ° " ^^^^ 

For a 7^ the expression is not more complicated but it is slightly uglier. A simple way to 
compute the path ordered exponential is to convert it into the solution to a linear problem. 
That linear problem is the weak coupling counterpart of the strong coupling linear problem 
[29]. That is, the path ordered exponential in (77) can be obtained from the monodromy 
matrix L{a) through VL{u) = L [log(S')/2] with 

[d. - uXcak(n; a)] L[a] = , L[- \og{S)/2] = J J ^ . (79) 

The linear problem can be easily solved and leads to some matrix Q{u) which we computed. 
That is, we just computed the right hand side of (77). This is an open holonomy so it is 

^^Naively it is not obvious that we are allowed to integrate up to log(S')/2 in (77), since the entries 
in the weak coupling connection (78) explode as S' — >■ oo. They naively invalidate (76). However the 
matrix (78) is not large in the sense that the small solutions do not diverge when acted with (78) many 
times repeatedly. An easy way to see this is the following. There exists a discrete gauge transformation 
Av,cak{u', cr) — T' A'^Jl]^{u; a) = g{a) ^wcak(M; cr) g~^{<^ + 1/j) such that (see section 5.2 for discussion) 

■ new I \ ^ f -«-2Mtanhcr 2uC0thcr \ r^r^ , ■2\ r n lon\ 

^wcak(";-)-^(^ 2.tanha z - 2. coth a j + ^^^/^ ) ' ^"'^ " = ^^^^ 

where g is such that 5(cr)(l + .Awcak(o'))(?~^(cr) is diagonal. This connection does not diverge and therefore 
we can compute the path ordered exponential by solving the linear problem in this gauge. This is equivalent 
to solving the linear problem in the previous gauge (where the connection appears large) and applying the 
gauge transformation in the end. Since it is legitimate to integrate the solution to the linear problem up 
to log(S')/2 in the new gauge and transform the result back to the previous gauge, we can simply do the 
computation in the original gauge and the result (81) is valid. Of course we also confirmed (82) numerically. 
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not a gauge invariant quantity. One possibility to construct a gauge invariant quantity, as 
discussed in section 5, is to contract this holonomy with two small solutions as in (59). 
Another simple option, presented in section 5.3, is to consider the folded situation where the 
back of the bottom and top loops are the mirrors of their fronts. Then instead of (59) we can 
compute the trace of the closed holonomy. The two options are related by (66). That second 
option is more convenient at strong coupling since it will correspond to a folded string where 
the two folds are symmetrical. That is, we compute 

front back 



-^closed — Tr 



Tr {n{u) ■ n{-u)-^) (81) 



To obtain the continuum limit (77), we need the R-matrices to change slowly, that is for the 
insertions to be inserted densely, with j ^ 1. Since the exponent in (77) scales as j/u we 
also need m ~ j 1 to have a good scaling. In this limit we find 



^closed (w) = „ , , + 2 

F(u) 



Piu) 



, 2 J Piu) log 5 , , 

cosh ^ ^ > ^ 82 



u 



where 

P{u) — '^'^ ~ ^ — ^ i tanh — (83) 

The transfer matrix is the trace of the monodromy matrix or, equivalently, it is the sum of 
its two eigenvalues e*^^"-* and e"*^*^"^. The functions p{u) are the so called quasi-momenta 
introduced for closed strings in AdS^ in [34]. Next we take log(S') to be large. Then Tdosed 
is exponentially large and the quasi-momenta reads 

p{u)-2 ^ log S (84) 

This is the main result of this section. 

The quasi-momenta encodes the various charges of the state. For example, the energy 
can be obtained by expanding it at m = [34] 

P'H = -^-^^ + ^W (85) 



E= 1. . (86) 



From (84) we then derive 

yiog'^ 

Jcosh^(a/2) 

As a consistency check we can recalculate this energy directly using our Hamiltonian repre- 
sentation (13). By acting with the Hamiltonian, one can easily read off the one-loop energy, 

r cosh a + cosh 4 i 2 n 

E = 2(2/) / da J log ^ = / J ,^A l + 0{l/j)], (87) 

J 1 + cosh a Jcosh (q;/2) 

-|log(5) 

The first 2 comes from the fact that the Wilson line is folded, 2g^ arises from the normaliza- 
tion of the spin chain Hamiltonian and j appearing in the integral is the density of insertions. 
As we see, the direct computation (87) and the algebraic curve one (86) agrees perfectly. 
Next, we move to strong coupling. 
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6.2 Strong Coupling 



In this section we will analyze the string solution dual to the strong coupling limit of the 
square Wilson loop with insertions considered in the previous section. 

The four cusp solution with no scalar insertions is related by an analytic continuation to 
the GKP string [36] which is dual to a twist two single trace operator at large spin [37], see 
also appendix B of [1]. Similarly here, in the presence of the scalar insertions and at a = 0, 
that solution is related by an analytic continuation to the folded string solution dual to a 
single trace operator at large twist and spin. That solution was studied before in [37-39] (see 
also references therein) in the context of the spectrum. The solution at a 7^ is included in 
a family of solutions that were studied in [40] in the context of polygon Wilson loops. 

The solution can be embedded in an AdS^ x subspace. The embedding coordinates 
parametrizing the AdS and the circle obey 

X\ + - Xl -Xl = \ and + Y^ = 1 (88) 

correspondingly. In these coordinates the solution reads 

( X Vx ' ^X + X ') " e*"^™*'^e'^"2e""=' = g and Fi + 1Y2 = e""''^'""-" (89) 

The original Alday-Maldacena solution, without any R-charge, corresponds io 9 = = ti / A. 
Here we used the notation a for one of the worldsheet coordinates because it coincides with 
the a direction introduced above (73). This was one of the two symmetries of the reference 
square. The other symmetry was translations in the OPE time t. It almost corresponds to 
the other worldsheet coordinate, namely t = t cot if. The range of a is related to the spin 
S as a E [— I logs', ^logS*]. The frequency Wt is directly related to the angular momentum 
in the sphere which corresponds to the total R-charge which is flowing, namely J. More 
precisely, 

Finally, we can relate the parameters appearing in AdS and in the sphere through the 
Virasoro constraints. We find 

2 fi^ + f) 2cos(2g)cot((^) 

^ = ilFTwM ' = j • 

The quasi-momenta are the eigenvalues of the total monodromy, 

U-^ ■ n{x) ■ n-\-x) ■ U = diag (e'^'("\ e-'P^^^) (92) 

where [34] 

ilogS 



da Aa(x) 



A. = ^ . 93 

1 — 



^^Note that the embedding of the Euchdian solution (89) in the sphere is complex. It is related by analytic 
continuation to a real rotating timelike solution. 
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Here, the Zhukowsky variable x is related to the spectral parameter u as 



xiu) = ^ 94 

It is yet another convenient parameterization of the spectral parameter. As usual, the prod- 
uct of the two monodromies in (92) comes about due to the folded nature of the string. In 
practice they just yield an overall factor of 2 for p{x). As in the previous section, these 
holonomies can be easily computed and p{x) can be read off straightforwardly. Let us now 
discuss the results and, in particular, the comparison with weak coupling. 

First let us consider the simplest case where 6 = /A. In this case the string is point like 
in the sphere since Wa = 0, see (91). For this case, the quasimomentum reads 



2x -f 

x—ir^4 



p(x) = — \ll + '--x^ \og{S) for e = n/A. (95) 



There are two limits of this expression which we would like to comment on. The first is the 
large u limit where 

4:7TU 

x{u) ~ — > 1 (96) 

This limit corresponds to the Frolov-Tseytlin limit [13,33]. This is the relevant limit for 
large charge j ^ 1. In this limit we can drop the I's in the quasi-momenta, that is 

e = 7r/4 
x/') fixed 



p{x) ^ '^^ J ~ uHog{S) for <( j>l . (97) 



This matches precisely with the weak coupling result (84) for a = 0. The second interesting 
limit of (95) is when the charge j — )■ is small. In that limit 

p{x)^--^^^\og{S), for {-^l^^ ■ (98) 

This form matches precisely the behavior of the strong coupling Y-functions in the limit 
where the the conformal cross ration in the x~ direction is large. That is, under the iden- 
tification logxi342 = —2 log 5* and 9' = ^ log in (71) we get that 2 sinh(6'') log xr342 = 

Finally, let us move to the general case with generic 6. The expression for arbitrary x 
and j is not particularly illuminating^^ but it simplifies quite a lot in the Frolov-Tseytlin like 
limit where j and x are large (with fixed ratio). Then 



p(x) ~ ^ log S 99 

X 

Using (96), this matches precisely the weak coupling result (84) provided we identify 

cos(2^) = -tanh(a) . (100) 
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It is given by p{x) 



V2x log(S) csc{ip)^{x^ + l) cos{2ip)-x^~2ix cos(2e) sin(2ip) + l 
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While all these matches are great, we should of course be very clear and keep our original 
motivation in mind: the main point of all these sections is not the precise match of the weak 
and strong coupling holonomies. After all, we know that these kinds of matches are accidental 
and very particular to the near BPS nature of the Frolov-Tseytlin limit. Instead, the main 
motivation for studying the continuum limit was to make sure we had correctly identified the 
weak coupling analogue of the strong coupling Integrability structures. That claim merited 
further justification due to the striking difference between the weak (69) and strong (71) 
coupling results. We claimed that the very different looking structures were not unexpected 
since in one case we were dealing with a few flux tube excitations (weak coupling) whereas 
in the other case we had infinitely many particles (strong coupling). For justifying the claim, 
we showed that the weak and strong coupling structures become much more similar to each 
other once many excitations are flowing already at weak coupling. Checking this statement 
was the goal of these sections and indeed, it was clearly fully confirmed, to an even larger 
extent that we really needed! 



7 Discussion and Future Directions 

In this paper we studied the OPE for null polygon Wilson loops at weak coupling. We mapped 
the computation of the expectation value of the loop to a sum of two point functions of local 
operators. A lot is known about the Integrability properties of local operators and this 
allowed us to translate back this knowledge to the Wilson loops. For example, up to now, it 
was hard to predict the leading OPE discontinuities for examples where multiparticles states 
play a role. From the local operator point of view, the number of particles is simply related 
to the length of the corresponding spin chains. From this point of view, different lengths 
are equally easy to deal with and hence we can tame the multiparticle example in the OPE 
and predict leading OPE discontinuities for any amplitude. This was the subject of the first 
sections (2, 3 and 4). 

Of course, ultimately, we want to find an alternative, Integrable, approach towards com- 
puting the full amplitudes at any value of the coupling and not only their leading disconti- 
nuities. To progress in this direction we went to strong coupling for inspiration. There, the 
key objects are holonomies of the flat connections between edges of the null polygons. These 
holonomies are very natural in Integrable models; they are generating functions of all the 
conserved charges. They contain information about the energy of the state - which was the 
key ingredient for the OPE discontinuity story - but they also contain information about all 
the other higher charges. We identified the analogue of these objects at weak coupling and 
computed them for a few examples. This was the subject of the last sections (5 and 6). 

There are now two obvious pressing questions: 

• We have the holonomies at weak and strong coupling. How to compute them at higher 
loops at weak coupling'^° or, rather, at finite coupling? Can we make an educated 
guess for their finite coupling form? Note that there is a major difference between 

'^'^Note that higher loop R-matrices are not known. Can one construct them? We have a lot of data we 
can use from the studies of amplitudes. Can this data be used to learn about these higher loop objects? 
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perturbative computations and finite coupling. In perturbation theory we have, at each 
order, a maximum number of particles propagating in the fiux tube. On the contrary, 
at finite coupling, there are always infinitely many excitations fiowing. From this point 
of view, we expect the strong coupling holonomies to be much more representative of 
the full quantum solution. '^^ 

• Note that computing the holonomies is not the end of the story. After all, what we 
are interested in is the expectation value of the Wilson loops! At strong coupling, the 
expectation value of the Wilson loop is the area of a minimal surface. This area is 
encoded in the asymptotic values of the holonomies. What is the analogue of this last 
step at weak coupling or at finite coupling?'^^ As mentioned in the previous point, 
it might be that the perturbative regime is too degenerate and it might be somehow 
simpler to guess the full quantum result in one go. 

We should have these big questions in mind but they are probably too hard and vague 
to tackle as stated. Hence, we end with some more pragmatic next steps. At weak coupling: 

• It would be interesting to consider more examples at weak coupling to gain some more 
experience and intuition. For example, using the techniques explained in sections 2- 
4, we can now compute the leading OPE discontinuities of N^MHV amplitudes. It 
might be interesting to bootstrap them following what was done for MHV and NMHV 
amplitudes. The main difference will be that for N^MHV amplitudes one is dealing 
with two-particle states while for MHV and NMHV single particle states were enough. 
We already computed one example of an OPE discontinuity for one N^MHV example 
at one and two loops, see (16) and (43). 

• The two main ingredients in the OPE are the energies of the fiux tube states and the 
probability amplitudes for creating and absorbing these states. The latter are denoted 
as the OPE form factors. As discussed in section 3, they receive contributions from the 
geometrical expansion of the top and bottom parts of the loop but also from quantum 
corrections due to interactions between the insertions. Can the latter corrections be 
taken into account by a spin chain operator like the spin chain Hamiltonian? In princi- 
ple, at leading order in the couphng, it would be straightforward (and very interesting) 
to compute these quantum corrections. Technically, this should be very similar to the 
computations in [45] and [46]. 

■^^At strong coupling, the holonomies often enter with shifts in the spectral parameter. These shifts are 
known as crossing shifts under which the Zhukowsky parameter x(u) (94) transforms as x ^ 1/x. That 
transformation does not commute with perturbation theory (see [41] and [3, 42] for its study in the context 
of polygon Wilson loops). Therefore, one may need to go to finite coupling, where crossing is not degenerate, 
in order to reveal the all loop structure. The anomalous dimension of a Wilson loop with a cusp was recently 
studied in [43]. In particular, in [44] an effective description of the full quantum system emerged where the 
fundamental shift in the spectral parameter is the crossing shift. It suggests that there may by a description 
of polygon Wilson loops where the crossing shift is the fundamental one at any value of the coupling. 

■^^Should we first build Y-functions out of ratios of holonomies? And, if so, how should we think of these 
ratios and products once they involve holonomies that cross? At the operatorial level with quantum R- 
matrices at the crossing points or at the level of products of classical expectation values? We suspect the 
former. 
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• So far, the power of supersymmetry has not been combined with the OPE (apart from 
a simple use of SUSY Ward identities in [6]). That is, we normally consider particular 
components of the super Wilson loop. Supersymmetry relates the different components 
and therefore one should be able to package together different OPE results in a compact, 
manifestly supersymmetric form. The use of supersymmetry was very important in 
fixing the key ingredients of the spectrum problem solution, namely the dispersion 
relation and the world-sheet S-matrix. It might be that it will be equally important in 
the OPE. Related to this point, a remarkable alternative approach towards computing 
the all loop S-matrix in A/^ = 4 SYM was put forward recently by Caron-Huot and 
Song He [10] and by BuUimore and Skinner [11]. In this approach a recursive relation 
for the Wilson loop was derived based on the Yangian symmetry of the problem and 
supersymmetry plays an absolutely central role in this construction. Another instance 
where supersymmetry and the Yangian are manifest is in the Grassmanian formulation 
of scattering amplitudes [47]. It would be wonderful if one could establish further 
connections between the OPE and these two exciting developments. 

At strong coupling: 

• As mentioned above, we expect the strong coupling result to be a very good repre- 
sentative of the full quantum solution. The strong coupling result takes the form of a 
thermodynamic Bethe ansatz free energy. We interpret this free energy as describing a 
sum over densities of particles flowing in the OPE flux tubes in the different channels. 
In other words, the free energy is the full re-summation of the OPE. However, this is 
just an interpretation of the final result, the derivation is a purely geometrical one and 
did not rely on this physical picture in any way whatsoever! It would be extremely 
instructive to re-derive the strong coupling result from the OPE point of view. We 
expect such a derivation to have a rather direct finite coupling generalization. 

• In the spectrum problem the R-charge of the local operators was directly related to 
the length of the corresponding spin chains. Solving the problem for large lengths was 
the first step towards the full solution. At the end, one can take the zero R-charge 
limit and also study purely gluonic operators. "^'^ Furthermore, in a classical large length 
limit it is possible to make direct connections between weak and strong coupling. In 
the classical limit both are described by algebraic curves [33] that become identical in 
the so called Frolov-Tseytlin limit [13]. Having this common language (the algebraic 
curves) was instrumental in guessing the quantum Bethe equations [48]. It is natural 
to expect that by adding large amounts of R-charge to Wilson loop, we might obtain 
a less degenerate situation where the weak and strong coupling descriptions become 
considerably more similar. This is exactly what we observed in section 6. We saw that 
the weak and strong coupling holonomies for a four-cusped polygon with R-charge 
insertion match precisely in an analogous limit. Hence, it might be simpler to first 
solve this problem by artificially adding a lot of R-charge that is then removed at the 
very end. As a first step, one could try to generalize the strong coupling Y-system to 
the case where we also have movement in the sphere. '^^ 

^^This idea was recently used in the quark-anti-quark potential problem in [43]. 

'''''In the Euclidean case, the sphere is complexified and it might be possible to think of it as another AdS. 



39 



We hope to come back to these points in the future. 
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A NMHV Octagon from the Super OPE in Momen- 
tum Space 

In the usual OPE approach for super loops [6] we start with an NMHV tree level amplitude 
and OPE promote it to higher loops using the dispersion relation 7(p) of the relevant particle 
propagating from the bottom to the top of the polygon. In the octagon Wilson Loop studied 
in section 4.1, see figure 10, the particle is a scalar and [3,4,49] 

7(p) = 2/ [^(1/2 + ipl2) + ^(1/2 - tp/2) - 2^(1)] (101) 

where %l){x) = r'{x)/r{x). Tree level amplitudes are most conveniently written using dual 
momentum twistors. In M^'^ these twistors can be parametrized as 

- \ T / + \ T 



Z2n-l OC [0, 0, 1, YZ^J ' « (^1, 0, Oj (102) 

Conformal transformations act linearly as SL(4) transformation of the twistors. We can act 
with the conformal transformation'^^ 



M 



[ T-y^ \ 



5+1/2 

V 5-1/2 J 



(103) 



on the bottom twistors. These generates a family of twistors where the bottom coordinates 
are changed according to x+ — )■ x^(T) and x~ — )■ x~{S) where xf{T) is defined in (6) and 
x~ {S) is defined similarly. 



The Y-system dispersion relation should now be modified. It should be something like in footnote 29. The 
two Y-systems should be coupled. 

■^^ Obviously, in this appendix S = e^°' parametrizes the OPE translations, one of the symmetries of the 
reference square. It should not be confused with the spin S used in most of the main text. There S was a 
cut-off, the total spin of the local operators. 



40 



To summarize, we can start with some random numerical twistors Zj. For the octagon 
we will have eight such twistors. Then we act with the transformation (103) on the bottom 
twistors. This generates a family of polygons W{S,T). There are only two independent 
cross ratios for the octagon so the two parameters 5* and T are all we need to describe all 
conformally inequivalent octagons. 

With this prescription all brackets (ijkl) which involve a bottom twistor will get T or S 
dependence. We have 

^(2367) ^ I ^ i . 

'^tree ^2367) 5-1/2 + ^+1/2' ^"^^^^ 

for a particular choice of the initial twistors. '^"^ Other choices of the initial twistors would 
only differ by translations of log(T) or log(S'). Next we Fourier transform this result with 
respect to a = — ^ log S to get 

The OPE promotion is now trivial. We simply have to include powers of the anomalous 
dimension (101) to get [1,6] 

2^^^ . J .(P) ^ 2/^<r' .0. (106, 



^ 4^(2367) ^(2) _ f ^ 2. ^ _ 4 4^( 



which coincide precisely with the predictions (36) and (40) derived in the main text using 
the Hamiltonian insertion method! 



B Hamiltonian Density Averages 

In this appendix we quote some results for Hamiltonian averages {'Ha+i'Hjj+i) . As discussed 
in section 4.1.2, there are two type of averages which are non-trivial, see figures 11 and 12. 
We will present the results for the two cases only in the kinematic regimes needed for the 
two loop examples discussed in the main text. To compute any of these examples we apply 
(11) once followed by (13), see main text for more details. For we find 

C • I 



{Zg^z^jn'oZa^z,)) ^ @ ^ (x+-i)Li2x+-(x'-i)Li2x" 



[Zd Zc Za-^ Zh) 



a I \ 1 



(107) 



^^The twistors are now chosen to be Zi = (0,0, 1,0)^ Z2 = M(l, 0, 0, O)-' , Z3 = M(0,0,1,1)^ Z4 = 
M(l, 1,0,0)^, Z5 = (0,0,0,1)^, Ze = (0,1,0,0)^, Z7 = (0,0,1,-1)^ and = (1,-1,0,0)'^. With this 
parameterization, the cross ratios are xr243 ~ ^ 5' X1243 = ^ + T. 
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where is given in (14). It is quite remarkable that the difference from the one loop average 
(13) is a simple replacement of a logarithm by a dilogarithm! Probably a similar replacement 
will give us the higher loop results 

The averages of {'Hi+2,i+i'Hi,i+i) will depend on the six cross-ratios (3 in the + direction 
and 3 in the — direction) constructed out of the three points in the bottom and the three 
points on the top which these Hamiltonians probe. There is no conceptual difficulty in 
computing ("Hj j+i'Hj+i 4+2) in general M}'^ kinematics. Instead, we will only present the 
cases needed for examples in the main text. For these cases some of the points are null 
separated. 

For the Octagon N^MHV example (43), the ("H^^i+i) result (107) reduce to 

*d c, 

C • ( 1 ■ J 

SI 1 / " 



- ' • - -2Li2(l - x,f,J , ^ = -2Li2(l-xIfeJ 



[Zd Zc Zg-k Zb) " ' 



(108) 

where the cross ratios are defined as in (14). The other cases need for the example (43) are 

d /I/ = ^^2(1 - Xabde) ,d /I 

/ ^ V /' 



log Xacef log Xabde 



/ \ • 

a ! \h 



\ I 



n) ! =l0gXa6del0gXfece/ \ = log Xabde log Xfcce/ 



/ n 



+Li2(l-X6cJ (h) \ +Li2(l-X 



^abefJ 



lb \ c / \b ' c 



(109) 

The notation should be self-explanatory: a gray line uniting two points indicates that they 
are null separated (not necessarily adjacent). So for example in the second figure we have 

For computing D^^ in (43) we need (top|EI^|bot)MHV5 (top|EI|bot)MHV, (top|EI|bot) and 
(toplM^lbot) = (top| (7^4,1 + "^1,2 + ^2,3 + ^3,4) ^i,2|bot) + cyclic (110) 
Using (108) and (109), we obtain (43). 
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C Non-compact Spin Chain Hamiltonians 



As explained in the main text, in the OPE approach the top and bottom parts of null 
polygons can be identified with the single trace operators made out of covariant derivatives 
and all other fields. 

The covariant derivatives are always present. They create the extended Wilson lines. We 
can also have other fields such as the scalars and the fermions in TV = 4. Which fields are 
involved depends on which excitations are being considered in the super Wilson Loop. In 
other words, it depends on the helicities of the particles being scattered in the dual Scattering 
Amplitude picture. 

The single trace operators can be thought of as non-compact spin chains. For the exam- 
ples considered in this paper two kinds of single trace operators appear: SL{2) operators 



rTr 



ni! 



D1^Z...DTZ 



and SL{2) x SL{2) operators 

|(ni,mi), . . . , {nL.niL)) = 



1 



riilmil . . . ni}.mi}. 



rTr 



(111) 



(112) 



where Z is a complex scalar. The integers ni,mj = 0, . . . , oo so we are dealing with non- 
compact spin chains where at each site we have an infinite number of states. 

One small parenthesis on the scalar products of two spin chains: Consider a ket \ni, . . . ,ni) 
representing an operator at (1,1) and a bra (mi, . . . , m^l representing an operator at (x"*", x~) = 
(0,0).'^'' Their scalar product is inherited from the tree level Wick contractions. More pre- 
cisely we write (mi, . . . , m^jni, . . . , n^) to indicate that the ath field of the first single trace 
operator is Wick contracted with the ath field of the second operator. Then we should take 
ni + mi derivatives of the first propagator - leading to a (— l)™^(?2i + mi)! factor - then 
n2 + 1712 derivatives acting on the second wick contractions and so on. Hence'^* 



(mi, . . . , m^lni, . . . , ni) = JJ(-l)'"" 



{na + niay. 



a=l 



na\m, 



a- 



It is sometimes useful to preform a change of basis where the states are orthonormal, see 
footnote 12. 

To OPE promote Null Polygon Wilson loops we need to know how to act on these states 
with the A/" = 4 planar dilatation operator. By planarity, the one loop dilatation operator D 
is a sum of Hamiltonian densities Hj^j+i which only act on the two sites i and i + 1, 



(113) 



■^^Recall that in the approach outhned in the main text we always expand the bottom and top Wilson 
lines around the points (0,0) and (1, 1) respectively. 

■^^We are using {Z{x' ,y')Z{x,y)) ~ . Of course there could be a constant multiplier but in all 

the physical quantities that we will consider it would drop out. 
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The two loop dilatation operator involves interactions between the three sites i,i + l,i + 2 and 
so on. For the purpose of the current paper all we need is the one loop dilatation operator. 
That is, we need the action of the Hamiltonian densities on SL{2) and SL{2) x SL{2) states, 

m\n,n') and H|(n, m), (n, m')) . (114) 

This action can be read off from the works of Beisert [23]. In fact, from these papers one 
can extract the action of the one loop dilation operator on any single trace operator. 

In this appendix we review what the action (114) is and derive some useful representations 
that turn out to be very convenient for our OPE purposes. In C.l we discuss the action on 
SL{2) states while in C.2 we consider the more involved action on SL{2) x SL{2) states. 



C.l Action on SL{2) states 

Here we discuss the action EI|n, n'). We will first present the action on states where the first 
site is empty. As we will review below, this is enough to fix the full Hamiltonian. 

Action on states with nothing on the first site 

The action of the one loop Dilatation operator can be neatly written using the harmonic 
numbers as 

n " 1 

m\0,n) = h{n)\0,n) -'^-\d,n-d) , /i(n) = ^ - (115) 

d=i j=i ^ 

and h{0) = 0. There is a very elegant and equivalent way of packaging the relations (115) 
using a generating function 

W[x\y] = Z{x) e(^-^)^+Z(x) . (116) 
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as 

r 1 J J 

mW[0\y] = I -^—^(^W[0\y]-tW[0\ty]-{l-t)W[ty\y]) (118) 

In fact, in practice it is often simpler to derive (118) by a direct Wilson loop computation 
and from that read off the action (115) on local operators. 

Action on generic states 

From (115) or (118) we can easily derive the action on a generic state \n,n'). Indeed, the 
Hamiltonian commutes with the total operator 

D+: D+\n,n') = {n + l)\n + l,n') + (n' + l)\n,n' + 1) . (119) 

■^^When checking or deriving this relation is is useful to know about a simple integral representation for 
the Harmonic numbers which can be derived as 

Hn) = E - = E / = / dt^^ (117) 
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Hence, we have the recursion relation 

m\n + 1, n) = (^D+m\n, n) - {n + l)e|n, n + 1)) . (120) 

Combined with the initial condition (115) this yields the SL{2) Hamiltonian [23] 

m\n, n') = [h{n) + h{n')] \n, n') - ^^\n- d, n' + d) (121) 

where the sum goes over d = —n', . . . ,n except d = 0. An even simpler way of deriving this 
result is by using the action on Wilson lines (118). Translating the full expression by x (and 
relabeling y + x ^ y) we have 

mW[x\y]= -(w[x\y]-tW[x\x + t{y-x)]- (l-t)W[x + t(y-x)\y]) (122) 

Jo — 1) V / 

Expanding each Wilson line 

oo CX) 
n=0 n'=0 

and collecting powers of x and y in each side of (122) directly leads to (121). 
C.2 Action on SL{2) x SL{2) states 

Next we move to the more complicated and considerably less studied case of EI|(ri, m), {n', m')). 
As before we start by ]HI|(0, 0), (n, m)) which is actually all we will use in the main text. This 
sector is not closed, that is the action of the Hamiltonian density on SL(2) x SL{2) states 
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is not a linear combination of SL{2) x SL(2) states (112) alone.^^ 

m Z [DID^'Z) = h{n + m)Z [DID'^Z] (124) 

^^Sn'+m'>ol ,j f A B{n' + m',l + n + m-n' -m') X 

n'=Om'=0 V / V / 



xiDlD"^ Z] {DT^'D"'-"' Z 



m min(n— n',m— m') 



n'=Om'=0 1=1 



~J2Y1 ^n'+m'>0 I - 1 I - 1 I - - 1 I - - 1 X 



n \ j m \ j m W n 
n'j \m'j \m' + l)\n' + 1 



X B{n' + m' + 1,1 + n + m - n' - m' - I) [d^D"!^' D[Z^ [dI-""' D"' DI^Z 



n m mm{n-n' ,m-m') / \ / \ / \ / 



n'=Om'=0 i=l 



n' J \m' J \m' + / / \n' + 
X B{n' + m' + /, 1 + n + m - n' - m' - /) [d^D'!!' D[Z^ j^^n-n'-i^m-m'-z^z 

where -B(a, 6) = r(a)r(6)/r(a + 6) and Dz and -D^ are the other two null direction. Schemat- 
ically, 

D+ = Do-^i, D- = Do + i?i, Dz = D2 + tDs, D, = D2-tDs. (125) 

The simplest way to obtain (124) is to read it from the Harmonic action of [23].'^^ The 
expression (124) looks scary and hard to work with. Fortunately, we will only need the first 
two lines of this expression which we will furthermore bring to a more useful form for our 
OPE applications. 

C.2.1 Projected Hamiltonian action 

First note that we consider polygons in M}'^. The fields in the bottom and top part of the 
polygon are therefore displaced in the x+ and x_ directions but not in the or Xz directions. 



'"^We can use the equations of motion to relate DzD^Z to D+D^Z which means that we can choose a base 
of operators where or Dz never appear in pairs acting on the same site (if they would we could replace 
them by )'^"™'"''' °f p'^"''. With this choice, we can have at most Uj's or D^'s in a given site; never 

both. With this prescription the translation between the oscillator language of [23] and covariant derivatives 
is unambiguous. 

"^^In this language one first maps the covariant derivatives to bilinears of oscillators. Then the Hamiltonian 
can hop these oscillators from the first site to the second site and vice versa. When there is no hopping 
the corresponding weight is an harmonic number of (half) the total number of oscillators N — 2n + 2m. 
That gives the first line in (124). The other lines in (124) correspond to some (71-21) oscillators jumping from 
the second site to the first site. There are no jumps in the opposite direction (ni2 — 0) since the first site 
was empty to start with. The corresponding weight is —B{n2i/2, 1 + N/2 — n2i/2), see (F.5) in the second 
reference in [23]. This explains all the B functions in (124). Finally, the binomial coefficients that multiply 
these B functions come from the number of ways of choosing which oscillators jump from the second to the 
first site. 
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Hence, suppose we act with a single Hamiltonian on the top. That generates fields which 
contain derivatives in the z and z directions, the last terms in (124). When contracting with 
the bottom we will end up with propagators like 

[d\D\D"_^z) [x+,a;_,x„x,] [oiD^ z) y_, y.-]) (126) 

where we used x = x^^"^^ and y = x^^°^^°"^\ To leading order in the coupling this is simply 
proportional to 



dx+J \dx+ J \dy+ J \dy J ydx"^ J {x+ - y+){x^ - y^) - {x^ - y^){x-^ - y^) 

After taking the derivatives we end up with positive powers of {xz — Vz) in the numerator. 
Since at the end we should set Xz = Xz = yz = yz = 0, we get zero. Hence we can drop the 
last two lines in (124) if we are only interested in leading perturbative computations. 

At higher loops this argument no longer works since we can act on the top and on the 
bottom and easily end up with (dropping the + and — derivatives that play no role in these 
arguments) 

d d V 1 



dx"^ dyy (x+ - y+)(a;_ - y_) - (x^ - y^)(xg - yz) 

In this case the z and z derivatives are balanced and we no longer get zero. Sometimes, even 
for some higher loop examples, the polygon kinematics are such that we can still drop the 
last terms in (124). 

Given this discussion, from now on we drop the last terms in (124) and consider the 
projection of the Hamiltonian action into the SL{2) x SL{2) sector. 



(^:^"'^) l5L(2)x5L(2) = + m)Z (DID-Z) (127) 
-EE L' U' B{n + mA + n + m-n'-r^){DiD-'z)[Dr-'D-~-'z). 

n'=lm'=l V / V / 

C.2.2 Wilson Line Action 

In this section we re- write (127) as a Kernel action on a Wilson line connecting two points 
in M^'\ 

W[x', y'\x, y] = Z{x\ y') e^—')D-+{y-y')D+ y') . (128) 
That is, we want to generalize (118) to this more involved case and compute 



mW[0,0\x,y] = y ^^e(ZD!^D™Z) (129) 

nl ml 



n.m=0 



using (127). Since the computation is fun, we will do it in detail; the final result is (133). 
When plugging (127) in (129) we have two terms coming from the first and second lines in 
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(127). The first one is trivial to deal with following the same steps as in SL{2), section C.l. 
Hence we consider the contribution from the second line in (127): 

„„ n m / \ 2 / \2 
X y / \ 



-E 



nl m 

n,m=0 n—Om—O 



n rn / \ z / \ z 

E E C) iZ) B{n'+m',l+n+m-n'-m')[DtD-'z) {d^-' D^^-' Z 

r,l — t\ml—n \ / \ / 



Note that we dropped the constraint 5„/+m'>o which imposes that there is some hopping. 
That means that at the end we should take care of subtracting the contribution of this 
expression with n' = m' = which is quite simple to do. Let us now simplify this expression. 
First we use 



B{a,b) 



dtf 



a-li 



and introduce n = p + n' and m = g + m' to get 



dt X'" y 

p,q,n' ,m'=0 



n' 



m' 



t) 



p+g 



(130) 

(D^'D^'Z) (DlDlZ) 



n'\m'\p\q\ 



Now we use that 



a + h 
b 



dz (1 + zY+^ 
27ri z^+i 

where the contour goes around the origin and a, b are positive integers. We end up with a 
totally factorized integrand/summand 



^ dt 

T 



dz 



dw 



2'Kiz J 2'Kiw 



-X 



(131) 



where 
X = 



[x(l + z)tj [i/(l + w;)tj [x(l + -)(l-t)J [i/(l + -)(1 -t) X (132) 

p,q,n' ,m' =0 



W 



X 



(L)!^'D-'Z) (DlD^Z) 



n'lm'lplql 



= W[x{l + z)t,y{l + w)t][x{l + l/z){l - t), y{l + l/w){l - t)] 

Now we just need to take care of the remaining contributions. They are the n' = m' = 
term that we dropped above plus the first line in (127). They are quite trivial to work out 
to let us simply quote the final result 

dt 



mw[o,o\x,y] 



SL{2)xSL{2) 



t 



dz 



dw 



2TTiz T 27Tiw 



(133) 



X {W[0, 0|x, y] - W[x{l + z)t, y{l + w)t\x{l + l/z){l - t),y{l + l/w){l - t)]) 

We could now get the general case for the action on a state containing also derivative in 
the first site by either translating this expression or by using a recursion relation together 
with the initial condition (127). See section C.l for a more careful discussion in the SL{2) 
context. For our purpose (127) or (133) is all we need since we can always translate a pair 
of points so that the origin coincides with one of them. 



48 



References 



[1] L. F. Alday, D. Gaiotto, J. Maldacena, A. Sever and P. Vicira, "An Operator Product 
Expansion for Polygonal null Wilson Loops," [arXiv: 1006. 2788 [liep-th]]. 

[2] L. F. Alday and J. M. Maldacena, "Comments on operators with large spin," JHEP 
0711, 019 (2007) [arXiv:0708.0672 [hep-th]]. 

[3] B. Basso, "Exciting the GKP string at any couphng," arXiv: 1010.5237 [hep-th]. 

[4] D. Gaiotto, J. Maldacena, A. Sever and P. Vieira, "Bootstrapping Null Polygon Wilson 
Loops," JHEP 1103, 092 (2011) [arXiv:1010.5009 [hep-th]]. 

[5] D. Gaiotto, J. Maldacena, A. Sever and P. Vieira, "Pulling the straps of polygons," 
arXiv: 1102.0062 [hep-th]. 

[6] A. Sever, P. Vieira and T. Wang, "OPE for Super Loops," JHEP 1111, 051 (2011) 
[arXiv: 1108. 1575 [hep-th]]. 

[7] L. J. Dixon, J. M. Drummond and J. M. Henn, "Analytic result for the two-loop six- 
point NMHV amplitude in N=4 super Yang-Mills theory," JHEP 1201 (2012) 024 
[arXiv:1111.1704 [hep-th]]. 

[8] L. J. Dixon, J. M. Drummond and J. M. Henn, "Bootstrapping the three-loop hexagon," 
JHEP 1111 (2011) 023 [arXiv: 1108.4461 [hep-th]]. 

[9] P. Heslop and V. V. Khoze, "Wilson Loops @ 3-Loops in Special Kinematics," JHEP 
1111 (2011) 152 [arXiv:1109.0058 [hep-th]]. 

[10] S. Caron-Huot and S. He, "Jumpstarting the All-Loop S-Matrix of Planar N=4 Super 
Yang-Mills," arXiv:1112.1060 [hep-th]. 

[11] M. BuUimore and D. Skinner, "Descent Equations for Superamphtudes," 
arXiv:1112.1056 [hep-th]. 

[12] L. F. Alday, J. Maldacena, A. Sever and P. Vieira, "Y-system for Scattering Amplitudes," 
J. Phys. A A 43, 485401 (2010) [arXiv:1002.2459 [hep-th]]. 

[13] S. Frolov and A. A. Tseytlin, "Rotating string solutions: AdS / CFT duality in non- 
supersymmetric sectors," Phys. Lett. B 570, 96 (2003) [hep-th/0306143]. 

[14] A. V. Belitsky, "OPE for null Wilson loops and open spin chains," Phys. Lett. B 709, 
280 (2012) [arXiv:1110.1063 [hep-th]]. 

[15] L. F. Alday, B. Eden, G. P. Korchemsky, J. Maldacena and E. Sokatchev, "From cor- 
relation functions to Wilson loops," arXiv: 1007.3243 [hep-th]. • B. Eden, G. P. Ko- 
rchemsky and E. Sokatchev, "From correlation functions to scattering amplitudes," 
arXiv: 1007.3246 [hep-th]. • B. Eden, G. P. Korchemsky and E. Sokatchev, "More on 
the duality correlators/amplitudes," arXiv:1009.2488 [hep-th]. 



49 



[16] M. BuUimore and D. Skinner, "Holomorphic Linking, Loop Equations and Scattering 
Amplitudes in Twistor Space," arXiv:1101.1329 [hep-th]. • A. V. Belitsky, G. P. Ko- 
rchemsky and E. Sokatchev, "Are scattering amplitudes dual to super Wilson loops?," 
arXiv:1103.3008 [hep-th]. 

[17] L. F. Alday, J. M. Maldacena, "Gluon scattering amplitudes at strong coupling," JHEP 
0706, 064 (2007). [arXiv: 0705. 0303 [hep-th]]. 

[18] G. P. Korchemsky, J. M. Drummond, E. Sokatchev, "Conformal properties of four- 
gluon planar amplitudes and Wilson loops," Nucl. Phys. B795, 385-408 (2008). 
[arXiv: 0707. 0243 [hep-th]]. • A. Brandhuber, P. Heslop, G. Travaghni, "MHV am- 
phtudes in N=4 super Yang-Mills and Wilson loops," Nucl. Phys. B794, 231-243 
(2008). [arXiv:0707.1153 [hep-th]]. • Z. Bern, L. J. Dixon, D. A. Kosower, R. Roiban, 
M. Spradlin, C. Vergu and A. Volovich, "The Two-Loop Six-Gluon MHV Ampli- 
tude in Maximally Supersymmetric Yang-Mills Theory," Phys. Rev. D 78, 045007 
(2008) [arXiv:0803.1465 [hep-th]]. • J. M. Drummond, J. Henn, G. P. Korchemsky 
and E. Sokatchev, "Hexagon Wilson loop = six-gluon MHV amplitude," Nucl. Phys. B 
815 (2009) 142 [arXiv:0803.1466 [hep-th]]. • N. Berkovits, J. Maldacena, "Fermionic 
T-Duality, Dual Superconformal Symmetry, and the Amplitude/ Wilson Loop Connec- 
tion," JHEP 0809, 062 (2008). [arXiv:0807.3196 [hep-th]]. 

[19] L. J. Mason, D. Skinner, "The Complete Planar S-matrix of N=4 SYM as a Wilson 
Loop in Twistor Space," JHEP 1012, 018 (2010). [arXiv: 1009.2225 [hep-th]]. 

[20] S. Caron-Huot, "Notes on the scattering amphtude / Wilson loop duahty," 
[arXiv:1010.1167 [hep-th]]. 

[21] T. Adamo, M. BuUimore, L. Mason and D. Skinner, "A Proof of the Supersymmet- 
ric Correlation Function / Wilson Loop Correspondence," JHEP 1108 (2011) 076 
[arXiv:1103.4119 [hep-th]]. 

[22] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, "Dual superconformal 
symmetry of scattering amplitudes in N=4 super- Yang-Mills theory," Nucl. Phys. B 828 
(2010) 317 [arXiv:0807.1095 [hep-th]]. 

[23] N. Beisert, "The complete one loop dilatation operator of N=4 superYang-Mills theory," 
Nucl. Phys. B 676 (2004) 3 [hep-th/0307015]. • N. Beisert, "The Dilatation operator 
of N=4 super Yang-Mills theory and integrability," Phys. Rept. 405, 1 (2005) [hep- 
th/0407277]. 

[24] A. V. Belitsky and A. V. Radyushkin, "Unravehng hadron structure with generalized 
parton distributions," Phys. Rept. 418, 1 (2005) [hep-ph/0504030]. 

[25] H. Elvang, D. Z. Freedman and M. Kiermaier, "Solution to the Ward Identities for 
Superamplitudes," JHEP 1010 (2010) 103 [arXiv:0911.3169 [hep-th]]. 

[26] N. Beisert, "The su(2|3) dynamic spin chain," Nucl. Phys. B 682 (2004) 487 [hep- 
th/0310252]. 



50 



[27] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, "Generalized unitarity 
for N=4 super-amplitudes," arXiv:0808.0491 [hep-th]. • A. Brandhuber, P. Heslop and 
G. Travaglini, "One-Loop Amplitudes in N=4 Super Yang-Mills and Anomalous Dual 
Conformal Symmetry," JHEP 0908 (2009) 095 [arXiv:0905.4377 [hep-th]]. • H. Elvang, 
D. Z. Preedman and M. Kiermaier, "Dual conformal symmetry of 1-loop NMHV am- 
plitudes in N=4 SYM theory," JHEP 1003 (2010) 075 [arXiv:0905.4379 [hep-th]]. • 
N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo and J. Trnka, "Local Integrals for Planar 
Scattering Amplitudes," arXiv: 1012.6032 [hep-th]. 

[28] N. Beisert, C. Ahn, L. F. Alday. Z. Bajnok, J. M. Drummond, L. Freyhult, N. Gromov 
and R. A. Janik et al, "Review of AdS/CFT Integrabihty: An Overview," Lett. Math. 
Phys. 99 (2012) 3 [arXiv: 1012.3982 [hep-th]]. 

[29] L. F. Alday and J. Maldacena, "Null polygonal Wilson loops and minimal surfaces in 
Anti-de-Sitter space," JHEP 0911, 082 (2009) [arXiv:0904.0663 [hep-th]]. 

[30] N. Beisert, V. A. Kazakov, K. Sakai and K. Zarembo, "Complete spectrum of long 
operators in N=4 SYM at one loop," JHEP 0507, 030 (2005) [hep-th/0503200]. 

[31] N. Beisert, V. A. Kazakov, K. Sakai and K. Zarembo, "The Algebraic curve of clas- 
sical supcrstrings on AdS(5) x S**5," Commun. Math. Phys. 263, 659 (2006) [hep- 
th/0502226]. 

[32] L. D. Faddeev, "How algebraic Bethe ansatz works for integrable model," hep- 
th/9605187. 

[33] V. A. Kazakov, A. Marshakov, J. A. Minahan and K. Zarembo, "Classical/quantum 
integrabihty in AdS/CFT," JHEP 0405, 024 (2004) [hep-th/0402207]. 

[34] V. A. Kazakov and K. Zarembo, "Classical / quantum integrabihty in non-compact 
sector of AdS/CFT," JHEP 0410, 060 (2004) [hep-th/0410105]. 

[35] S. Schafer-Nameki, "Review of AdS/CFT Integrabihty, Chapter II. 4: The Spectral 
Curve," Lett. Math. Phys. 99 (2012) 169 [arXiv:1012.3989 [hep-th]]. 

[36] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "A Semiclassical limit of the gauge / 
string correspondence," Nucl. Phys. B 636 (2002) 99 [hep-th/0204051]. 

[37] R. Roiban and A. A. Tseytlin, "Spinning superstrings at two loops: Strong-coupling 
corrections to dimensions of large- twist SYM operators," Phys. Rev. D 77 (2008) 066006 
[arXiv:0712.2479 [hep-th]]. 

[38] L. Freyhult, A. Rej and M. Staudacher, "A Generalized Scaling Function for AdS/CFT," 
J. Stat. Mech. 0807 (2008) P07015 [arXiv:0712.2743 [hep-th]]. 

[39] N. Gromov, "Generalized Scahng Function at Strong Coupling," JHEP 0811 (2008) 085 
[arXiv:0805.4615 [hep-th]]. 

[40] H. Dorn, N. Drukker, G. Jorjadze and C. Kalousios, "Space-like minimal surfaces in 
AdS X S," JHEP 1004 (2010) 004 [arXiv:0912.3829 [hep-th]]. 



51 



[41] R. A. Janik, "The AdS(5) x S**5 superstring worldsheet S-matrix and crossing symme- 
try," Phys. Rev. D 73, 086006 (2006) [hep-th/0603038]. 

[42] B. Basso and A. V. Belitsky "Luescher formula for GKP string," Nucl. Phys. B 860, 1 
(2012) [arXiv:1108.0999 [hep-th]]. 

[43] D. Correa, J. Maldacena and A. Sever, "The quark anti-quark potential and the cusp 
anomalous dimension from a TBA equation," arXiv:1203.1913 [hep-th]. • N. Drukker, 
"Integrable Wilson loops," arXiv: 1203. 1617 [hep-th]. 

[44] N. Gromov and A. Sever, "Analytic Solution of Bremsstrahlung TBA," arXiv: 1207.5489 
[hep-th]. 

[45] L. F. Alday, J. R. David, E. Gava and K. S. Narain, "Structure constants of planar N 
= 4 Yang Mills at one loop," JHEP 0509, 070 (2005) [hep-th/0502186]. 

[46] J. Plefka and K. Wiegandt, "Three-Point Functions of Twist-Two Operators in N=4 
SYM at One Loop," arXiv: 1207.4784 [hep-th]. 

[47] N. Arkani-Hamed, F. Cachazo, C. Cheung and J. Kaplan, "A Duality For The S Matrix," 
JHEP 1003 (2010) 020 [arXiv:0907.5418 [hcp-th]]. • N. Arkani-Hamcd, J. L. Bourjaily, 
F. Cachazo, S. Caron-Huot and J. Trnka, "The All-Loop Integrand For Scattering Am- 
plitudes in Planar N=4 SYM," JHEP 1101 (2011) 041 [arXiv: 1008.2958 [hep-th]]. 

[48] G. Arutyunov, S. Frolov and M. Staudacher, "Bethe ansatz for quantum strings," 
JHEP 0410, 016 (2004) [hep-th/0406256]. • N. Beisert and M. Staudacher, "Long- 
range psu(2,2|4) Bethe Ansatze for gauge theory and strings," Nucl. Phys. B 727, 1 
(2005) [hep-th/0504190]. • N. Beisert, B. Eden and M. Staudacher, "Transcendentality 
and Crossing," J. Stat. Mech. 0701, P01021 (2007) [hep-th/0610251]. 

[49] A. V. Belitsky, A. S. Gorsky and G. P. Korchemsky, "Logarithmic scaling in gauge/string 
correspondence," Nucl. Phys. B 748 (2006) 24 [hcp-th/0601112]. • L. Frcyhult, A. Rej 
and M. Staudacher, "A Generalized Scaling Function for AdS/CFT," J. Stat. Mech. 
0807, P07015 (2008) [arXiv:0712.2743 [hep-th]]. 



52 



